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SUMMARY 


The Streamtube Curvature Program (STC) for prediction of the pressure 
distribution about a nacelle at transonic speeds was developed by the General 
Electric Conq)any under contract number NASl-10804. Experience with the use 
of the program has indicated the desirability of modifying the program to 
extend its capabilities and provide inq)roved diagnostic Information in the 
event of computational difficulties. This report is concerned with the im- 
provements which have been incorporated in the STC program to enhance both 
its conq)utational and diagnostic capabilities. In Part I, specific descrip- 
tions are given of the revisions to more reliably handle the jet stream - 
external flow interaction at trailing edges. Also described are the augmented 
boundary layer procedures and a variety of other program changes relating 
to program diagnostics and extended solution capabilities. Part II consists 
of the updated User's Manual, and includes information on the conq)uter program 
operation, usage, and logical structure. 
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SECTION 1.0 


INTRODUCTION 


Aircraft have been designed with the NASA-developed supercritical wing 
to fly at cruise Mach numbers approaching one. The need for low installed 
drag and high drag- divergence Mach number nacelle installations is necessary 
to the success of this design. Design techniques are required to evaluate 
these nacelles on an isolated basis and then on an installed or integrated 
basis . 

For this reason, NASA has initiated a program to provide design infor- 
mation for low drag, high drag- divergence Mach number isolated nacelles 
suitable for use with advanced high bypass ratio turbofan engines. One ele- 
ment of such a program is the development of a method to predict the inviscid 
pressure distribution and flow field about an arbitrary axisymmetric ducted 
body at transonic speeds. This prediction technique provides the means to 
conduct parametric studies so that the nacelle design criteria can be evaluated 
to select configurations for further experimental investigations. The predic- 
tion technique also provides guidance during wind tunnel testing to develop 
nacelle shapes which minimize drag within given design restraints. 

The objective of the development of this computer analysis was the pre- 
diction of flow fields about Isolated nacelles at transonic conditions . The 
solution technique has been further specified to give accurate results con- 
sistent with the requirement of relatively short computing time per input case 
as conqjared to that required for a time dependent finite difference method of 
solution. The method utilized to compute the flow field was developed under 
contract NASl-10804 and is termed the .Streamtube Curvature Relaxation tech- 
nique (Reference 1) . 

The Streamtube Curvature (STC) method of solving external flows has not 
been discussed significantly in the literature; however, the method is a very 
natural one. For example, engineers frequently rely on one- dimensional com- 
pressible flow relationships for a first-order solution to ducted flows. The 
STC approach is similar except that a number of confluent streamtubes, with 
slightly different properties, are added together to obtain the total flow 
in the channel. Each streamtube is handled in much the same way as is the 
one streamtube in the one- dimensional problem. In the limit, as the size of 
the individual streamtubes approaches zero, the STC method satisfies the 
inviscid equations of motion exactly. 

While this program is a basic tool for the calculation of the overall 
forces on nacelles, it is well understood that the effects of boundary layer 
friction drag and displacement of the inviscid flow mxist be included to 
obtain accurate performance predictions. The displacement of the inviscid 
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flow effectively changes the body shape, thus altering the inviscid pressure 
distribution. The summation of pressure-area forces taken over the body can 
be seriously in error when this displacement effect is not included. The 
coupling of a boundary layer solution with the inviscid STC analysis allows 
inclusion of displacement effects as well as a complete evaluation of all 
nacelle forces, including friction and the prediction of boimdary layer 
separation. 

The accumulation of running experience with the STC program has indicated 
the desirability of modifying the program to extend its capabilities, inprove 
its reliability, and provide enhanced diagnostic information in the event of 
computational difficulties. The methods of analysis and the basic numerical 
procedures for computerization of the analysis are described in Reference 1. 
Part 1 of this report is concerned with recent modifications which have been 
Incorporated in the STC program to provide the above cited improvements . 

The conputer program source deck for the updated version of STC, together 
with a user manual are available from COSMIC (CoiH)uter Software and Informa- 
tion Center), Burrows Hall, University of Georgia, Athens, Georgia, 30601. 

The program is written in CDC Fortran 2.3 source language, with exception of 
three subroutines in COMPASS 1.1 language. The computer program has been 
checked out for the CDC 6600 machine. 
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SECTION 2.0 


MODIFICATIONS TO THE STC PROGRAM 


In this section, the details of the iiqjroved procedures inqilemented in 
the STC program are presented. Information on the conqjuter program operation, 
usage, and structure is described in the Users Manual (Part II). Since the 
bulk of the work involved specific program modifications, liberal references 
are made to this section of the report. 


2.1 JET STREAM - EXTERNAL FLOW INTERACTION 

To Increase the accuracy of the STC calculation in the region of a 
trailing edge, the following in^jrovements have been incorporated into the 
calculation procedure. 

• When the total to static pressure ratios on both sides of a trail- 
ing edge are subcritlcal, it is necessary that the static pressure 
on the two sides of the trailing edge are matched. This is 
accoiqillshed by an iterative adjustment of the flow rates, either 
on one or both sides of the trailing edge. The flow adjustment 
iteration logic was modified so that the flow adjustment is now 
performed in a separate iteration loop, and a record of the 
iterative adjustments is printed. 

• When a supercritical total-to-static pressure ratio exists at the 
trailing edge, a Prandtl-Meyer expansion fan (or compression wave 
approximated by the Prandtl-Meyer formulas) is computed to define 
the local change in streamline angle and static pressure. 

• The wake thickness distribution has been modified so that the 
wake displacement thickness provides a smooth transition between 
the blunt body and the "far" downstream flow field. 

In the implementation of the new logic, three options have been provided for 
user selection. A description of these options will be provided in the 
discussion which follows. 

First, however, to illustrate the nature of the flow in the trailing 
edge region, a fine grid calculation has been made for a low speed case in 
which the stagnation pressures of the flows above and below the trailing edge 
have been chosen to be unequal. Results for two levels of grid refinement 
are shown .in Figures 1 through 4. As shown in Figures 3 and 4, a spike in 
pressure occurs due to the abrupt change in angle near the trailing edge. 
Although not quite predicted by the nvnnerical procedure, the static pressure 
on the side with the lower total pressure (in this case the upper side) 
theoretically attains the stagnation value. At the trailing edge point on 
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Figure 1. Calculation Grid for Detailed Sharp Trailing Edge Solution, 
7 Grid Refinements, 
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Figure 2. Enlarged Plot of Calculation Grid in Trailing Edge Region, 


9 Grid Refinements. 











the lower side, the static pressure will match that of the upper side in 
order to satisfy the "Kutta" condition. Consequently, because of the higher 
total pressure in the lower stream, the stagnation condition will not be 
obtained and the angle variation along the lower dividing streamline must be 
continuous as depicted in Figure 5. This detail, however, is so fine that 
it is not readily determined in the numerical solution. Figure 6 shows a 
plot of calculated angles. 

This example helps to illustrate several interesting facts about trail- 
ing edge flow, namely: 

• As is well known when the stagnation pressures are equal on both 
sides of the trailing edge, the trailing streamline bisects the 
trailing edge wedge and the flow stagnates at the tip on both sides 
of the trailing edge. 

• When the stagnation pressure on one side of the wedge is somewhat 
higher than on the other side, the general character of the flow 
does not change markedly. Pressure spikes still occur on both 
sides of the trailing edge, although the theoretical flow solution 
is truely "singular" only on the side with the lower total pressure. 

• The width of the spikes (i.e., the length of near-stagnated flow) 
is very narrow and would be even narrower for smaller wedge angles . 
The calculations shown were carried out for a 22.6® wedge angle. 

The first of the trailing edge options now available in the STC program 
allows the user to select either the numerically computed or theoretically 
exact value of trailing edge pressure. This option, controlled by the input 
value of PDUM(2) , applies only to sharp trailing edges and also applies only 
to the stream with the lower total pressure. The normal (or default) option 
is to print the theoretical (i.e., the exact) values of trailing edge tip 
pressure and to use these same values in the calculation of the axial pressure 
forces . 

Figure 8 illustrates an example of a STC output when the theoretical 
value of the stagnation pressure is printed. The calculation grid for this 
example is shown in Figure 7a. These results may be compared to the alternate 
option of printing the (approximate) numerically computed values. Figure 9. 

The theoretical stagnation point condition on the upper side of the trailing 
edge generates a sharp downward spike in Mach number and an upward spike in 
pressure. Figure 8. Because the grid is relatively coarse, the high pressure 
region is averaged over a significant portion of the body and the error in 
the global momentum balance is large (20%). The global error is the sum of 
the entering momentum flux, the Integrated axial boundary pressure forces, 
and the negative of the leaving momentum flux. These global errors are 
normalized by the axial pressure force along the respective wedge surface and 
slip line boundaries. 

The "numerical" value of trailing edge pressure [PDUM(2)=0] shown in 
Figure 9 can be thought of as a representative average for the grid interval. 
Consequently, the computed global error is smaller. 


12 




otal Pressures. 




saaaSap ‘aiSuv 


Figure 6. Trailing Streamline Angle Distribution, 9 Grid Refinements, Detailed Sharp Trailing 
Edge Solution. 


Figure 7a. STC Grid for Jet Stream-External Flow Interaction Test Calculation 
Sharp Trailing Edge. 



Figure 7b. 


STC Grid for Jet Stream-External Flow Interaction Test Calculation 
Blunt Trailing Edge. 
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2.1.1 Prandtl-Meyer Expansion Fan Calculation 


In the STC algorithm a test is now made to determine if a Prandtl-Meyer 
expansion fan exists on either side of a trailing edge. For a convergent 
nozzle, this fan expands the flow from the pressure corresponding to unity 
Mach number to the pressure which is defined to exist on the other side of 
the trailing edge, Pg, shown in Figure 10. The theoretical value of Pg for 
a sharp wedge is the external stagnation pressure. Results obtained by 
using this option are shown in Figure 11. The alternate choice is to utilize 
the numerically calculated external pressure to define the Prandtl-Meyer 
expansion pressure, Pg, an exan^jle of which is shown in Figure 12. Except 
for the magnitude of the Mach number jump at the trailing edge, the solutions 
for the two options are not largely different. As stated earlier, this 
option applies only to sharp trailing edges. When the trailing edge is blunt. 
Figure 13, Pg is always taken as the numerically calculated external pressure. 


2.1.2 Wake Thickness Calculation 

The distribution of the wake thickness downstream of a blunt trailing 
edge has been changed from a linear to a quadratic form. The form used is: 


b = b + (b 


t. e. 


b ) 

oo' 



( 1 ) 


where b is the local wake thickness, x is the distance along the wake center- 
line from the trailing edge and boo, bj-^g^ and are defined in Figure 14. 

This form has the advantage that the length, can be chosen to provide a 

smooth continuation of the blunt body where, if the flow is subsonic, the wake 
edge streamline will be tangent to the effective body surface. Namely, at 
the trailing edge. 



tan a 


( 2 ) 


where a is the effective trailing edge wedge angle including the boundary 
layer displacement effects as shown in Figure 14. Differentiation of Equation 
(1) and substitution into Equation (2) gives: 


w tan a ■ 


(3) 


Note that if the wake centerline is curved, the streamline tangency condition 
at the trailing edge may still be maintained. Also, for small wedge angles, 
tana>0.1, the wake length is limited to twenty "effective" trailing edge' 
thicknesses. 
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PDUM (2) =0: P = numerically cal- 

culated TE pressure 

PDUM (2) =1: P = total pressure of 

® the external stream 



Figure 10. Options Available for Prandtl-Meyer Expansion 
Pressure, P^. 
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Downstream of the wake length, Si^, shown In Figure 14, a constant wake 
thickness, boo, is assumed. In the present procedure, b,, is taken as zero, 
but could be changed if desired. Far downstream of the trailing edge, it 
can be shown that the wake displacement thickness becomes equal to the momen- 
tum thickness. The momentum thickness is directly related to the friction 
drag and will not, in general, be zero. However, the effect of the finite 
versus zero downstream wake thickness on the body pressure distribution should 
be negligible for most cases. 

When a Prandtl-Meyer expansion fan exists on one side of the trailing 
edge, the effective flow convergence angle, a, downstream of the expansion 
fan becomes larger, as indicated in Figure 15, and consequently the wake 
length is reduced. However, if the user wishes the wake length to be based 
upon the wedge angle upstream of the expansion fan, this may be enforced by 
inputting PDUM(3)=1. 


2.1.3 Curvature at First Point Downstream of the Trailing Edge 

Curvature at supersonic points is evaluated by fitting a parabola to 
the point in question plus the two upstream points. However, at the first 
point downstream of a trailing edge, this procedure is modified. The three 
pieces of information used to define the parabola are the two points , the 
trailing edge point and the point in question, and the angle at the trailing 
edge just aft of the Prandtl-Meyer expansion. 

Again an option is provided to override this latter procedure, if 
desired, and to use the standard three point curvature formula even for points 
just downstream of the trailing edge. This is accomplished by setting 
PDUM(4)=0; an exanple result is shown in Figure 16. When compared to Figure 
13, some differences may be noted in the angle distribution along the lower 
slip line. However, the differences in the computed pressure and Mach numbers 
are not large for this case. Presumably, they would become more significant 
for a refined grid. 


2.1.4 Flow Adjustment Procedure 

One of the important factors in establishing an accurate flow solution 
is the proper determination of the relative flow rates in the several passages 
or, expressed differently, the determination of the relative split between the 
flow which passes over and under an immersed body. When the flow is subsonic 
at the trailing edge, this flow split is determined by the condition that 
pressure at the tip of the trailing edge must be a single value; i.e. , the 
two joining flows must have the same static pressure at that point. This rule 
may be applied to either sharp or blunt trailing edge cases (providing the 
flow is nearly stagnant behind the blunt edge) and to cases with finite 
boundary layer thicknesses. 


24 





Supersonic Jet Stream-External Flow Interaction. 
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As a word of caution, it should be noted that if the body shape is not 
well streamlined and a flow separation occurs upstream of the trailing edge, 
then the proper flow split can not be adequately determined. This is 
because a flow separation, which is not accounted for in the STC program, 
interacts with the inviscld flow and renders an accurate trailing edge 
static pressure prediction impossible. 

As an example of the procedure followed to obtain pressure closure at 
the trailing edge, the flow rate Iteration history for the Detailed Sharp 
Trailing Edge Solution (Figures 1 to 4) is shown in Figure 17. As illustrated, 
a flow iteration is performed for each grid refinement level. Errors which 
arise from an imbalance in trailing edge pressure are determined after the 
basic "inner” iteration tolerance is satisfied. The flow rate is then 
adjusted and additional inner Iterations are performed, again, tintil the 
inner iteration tolerance is satisfied. 

The trailing edge pressure closure error, plotted in Figure 17, is 
defined as fractional change in (the variable) flow rate necessary to achieve 
the same static pressure on both sides of the trailing edge. In this evalpa- 
tlon, the current distribution of streamline curvatures are used. (Of course, 
if the flow rate is changed, the streamline curvatures will then be somewhat 
in error, thus requiring additional inner iterations as mentioned above.) 

The user, through input specification, determines whether the channel flow 
rates above or below the trailing edge are to be varied to achieve the trail- 
ing edge pressure closure. If the flow rate in both the channels above and 
below the trailing edge are to be varied, the program will then hold the total 
flow constant and equal to the input value. (In this case, the flow rate 
shown in the "Kutta Iteration" printout is the flow rate below the trailing 
edge). 

For completeness, it should be pointed out that the trailing edge static 
pressure used in the flow adjustment logic is always the numerically con 5 >uted 
static pressure and not the theoretical stagnation value. This provides a 
consistent procedure for both sharp and blunt trailing edges, for both coarse 
and fine grid spacings, and for both equal and unequal total pressures of 
adjoining streams. (Notice that if the theoretical stagnation tip pressures 
were used for a sharp trailing edge in irrotational flow, then no pressure 
Imbalance could be determined ^ls a function of the flow split) . 

The above remarks apply to the case when the flow at the trailing edge 
is subsonic. If the flow is supersonic 'at the trailing edge, the static 
pressure equilibrium is obtained by an expansion fan or conqjression wave, 
whichever is appropriate. Furthermore, no pressure signal will be fed 
forward in this case to cause the flow rate to be changed. 

In the jet stream- external flow interaction calculations shown in 
Figures 11, 12, 13 and 16, the jet flow rate was set by user input rather 
than by the adjustment procedure which is internal to the program. The jet 
flow rate was selected, by trial and error, to satisfy the criteria that the 
Mach number on the internal boundary trailing edge point (just upstream of 
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Figure 17. Kutta or Trailing Edge Pressure Closure Iteration for the Detailed Sharp Trailing Edge 
Solution. The Right Graph is Enlarged to Illustrate Final Detailed Adjustments. 
Numbers Beside Points are the Grid Refinement Level and Inner Iteration Counter for 
which the Closure Error was Evaluated. 




the expansion fan) is unity. This was accomplished by changing the input 
value of the area, AO, which is used to compute the channel flow rate. The 
final selected value of AO was 2.99 corresponding to a flow coefficient of 
99.7%. (Thus, the selected input for these cases was: PT0=25, TTO=1050, 

A0=2.99, MACH0=1.0, VAKy=F. From these values and one- dimensional formulas, 
the channel flow rate computed by the STC program. The physical throat area 
was 3.0.) 

A computerized procedure for calculating the choked flow rate is also 
available. In this procedure the maximum flow rate is found Cfor the ’’assimed” 
set of curvatures) which will pass through the known area. This procedure 
is valid for the throat station of a convergent-divergent nozzle, where the 
wall is curved and the flow is partially subsonic and partically supersonic. 

In the straight wall convergent nozzle, the latter procedure is not valid and 
must be suppressed, as Indicated above, by setting VARY=F in the channel 
input data. Future additions to the STC program will be addressed to extend- 
ing the algorithm for use with straight walled convergent nozzles. 


2.1.5 Flow Field Convergence 

For the calculation grid Illustrated in Figure 7b typical iteration 
histories are shown in Figure 18. As indicated, the rate of convergence was 
found to be dependent upon the magnitude of the velocity jump in the slip 
line. For the highest jet pressure ratio of 5.2, convergence was not obtained. 

A large number of runs were made with different "correction equation 
deceleration factors." For example, RHOC and RHOW were varied from .66 to 
2.0. However, no significant improvement was obtained by the use of these 
factors. The use of the convergence factor CNVF was not tried. 

Another factor which was foimd to greatly affect the convergence rate 
was the axial spacing of the calculations stations. Subdivision of the grid 
to half of the spacing shown in Figure 7 led to solution divergence. 

The solution divergence was found to result from the omission of certain 
coefficient terms in the matrix equation for the streamline corrections. For 
flows with large vorticity, these terms become significant and, unfortunately, 
are not included in the present STC matrix solution procedure. 

2.2 BOUNDARY LAYER REVISIONS 


The basic boundary layer and separation calculation procedures defined 
in References 1 and 2 have been preserved intact in the current version of 
the program. Significant improvements, however, have been made in the imple- 
mentation philosophy and the associated program logic which supports the 
boundary layer and separation calculations. 
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Figure 18. Convergence Histories for Supersonic Jet with External Flow using 
- 8 Calculation Stations of the Grid shown in Figure 7b. 
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2.2.1 Wake Adjustment to Reflect Trailing Edge Boundary Layers 

In the original version of STC, wakes were allowed to exist in the 
regions downstream of blunt trailing edges. The base regions enclosed by the 
wakes were gradually reduced to zero as the solution proceeded downstream. 

In the updated version, entries in the STC wake table are now initialized at 
all trailing edges rather than just at those which have a finite thickness. 
This modification permits a redefinition of the wake table to allow displace- 
ment effects at both sharp and bltant trailing edges. 



Using the known surface geometry and the 6* and from the boundary layer 
calculation, the displaced coordinates and angles are defined as: 


Upper Surface 
^CU ^TEU ^ ds 


U 


^CU “ ^U ^U “^CU 


:( 4 ) 

(5) 


Z = Z - 6 sin <f) 

CU BU U ^CU 


( 6 ) 


Lower Surface 




^TEL ds 


(7) 
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( 8 ) 


^CL 

^CL “ ^BL ^CL 


AR^ - (Rcu “ 



The effective wake thickness at t! 

= \J 


C9) 

CIO) 

( 11 ) 

( 12 ) 

trailing edge is then: 

(13) 


2.2.2 Initiation of Boundary Layers on Collated Boundary Segments 

The STC program provides for definition of a boundary in discrete seg- 
ments, where continuity of the surface is maintained by matching the coordi- 
nates and angles at the junction points of the segments. The resulting 
surface is then collated under a single boxindary name internally for use by 
the program. For inlet configurations, with inner and outer contours joined 
at a stagnation point, the boundary .segments are collated under the boundary 
name of the inner surface. In all other collation situations, the boundary 
is combined under the name of the first upstream segment. 

Provision has been made in the STC program to allow the initiation of 
a boundary layer at the initial point of any collated boundary segment. For 
inlet geometries, this permits the user to run cases where the inner contour 
is "clean" and a boundary layer is calculated only on the external cowl 
surface. For boundary layers on surfaces which are segmented in the stream- 
wise direction, additional information must be supplied to reflect the up- 
stream history of the boundary layer (viz: 6*, 6 at the initial point). 

Within the STC program, the calculated boundary layer data are stored in a 
fixed length table area, along with other pertinent STC information (see 
Part II). When boundary layers are calculated on a large number of surfaces, 
the table area is often of insufficient size to permit storage of the boundary 
layer data. Using the above procedure to calculate boundary layers only on 
portions of given collated surfaces provides more effective use of the STC 
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table area. In this situation, care must be exercised when atten 5 )ting to 
factor displacement effects into the inviscid solution, since boundary layer 
information is not available upstream of the initial point. The resulting 
discontinuity in displacement thickness on the boundary could conceivably 
induce tmdesired pertubations of the STC inviscid solution. 


2.2.3 Separated 6* Calculation 

In the original version of STC, the boundary layer calculation was 
terminated if separation on a given surface was indicated. This condition 
caused several problems. Initially, the cumulative friction drag could not 
be calculated, since skin friction information downstream of the separation 
was not available. Finally, the boundary layer 6* downstream of the separa- 
tion point was assumed equal to the separated value. In general, this would 
tend to ignore the reattachment phenomenon and yield physically unreasonable 
values, particularly in both favorable pressure gradient zones and in the 
acceleration region downstream of a shock induced separation. 

In an attempt to partially alleviate these situations, the boundary 
layer calculation is continued downstream after flagging the separation. While 
the numerical magnitude of the displacement thickness may be somewhat in 
question, it is anticipated that these procedures may be more representative 
of the true physical situation. Upon restart, a separated flow warning is 
still Issued, as was done in the original inq>lementation of the STC program. 


2.2.4 Boundary Layer Initiation - Grid Refinement Level 

An input option has been provided to allow initiation of a boundary 
layer at any predetermined level of grid refinement. In many cases, it is 
desirable to obtain a fully converged inviscid solution before introduction 
of the boundary layer. The new input option allows a staged solution of 
the problem using the "restart" feature with the final tolerance reduced on 
the last restart. The boundary layer calculation may then be carried out 
at this point. With the standard form of input, boundary layers must be 
specified at the beginning of the problem, and the calculation will occur 
prior to the execution of the restart cases. The details for this input 
option are given in Part II. 


2.3 REVISED PROGRAM OUTPUT 


Several major program changes have been made in both the standard output 
and in the diagnostic output which occurs when the solution encounters com- 
putational problems. The ultimate objective of these modifications has been 
to provide the user with a clear and concise definition of the solution so 
that the accuracy of the results may be rapidly accessed and solution problem 
areas may be easily located. Detailed descriptions of the new output are 
given in Part II. A brief description is Included herein. 
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2.3.1 Iteration History 


The parameters In the Iteration hlstoiry output have been altered to 
provide a more meaningful picture of the solution development. The Items 
pertaining to the streamwlse point movement have been eliminated. New output 
Includes the coordinates of the point where the maximum flow balance error 
occurs. The algebraic sign of the maximum error has been retained so that 
the \iser may detect oscillation In sign as well as location In the field. 

The flow adjustment history at all trailing edges Is now printed. This 
Item Includes the X12 location, the current flow rate, and the flow adjust- 
ment error as related to satisfaction of the ’’Kutta" pressure con 5 )atablllty 
condition at the trailing edge. The final new Item of output Is the nxamber 
of Imbedded supersonic points In the flow field. This parameter provides a 
useful Indicator of what fraction of the field Is transonic and often points 
to possible solution problems. 


2.3.2 Diagnostic Format Revisions 

The format of the diagnostic output, produced when solution problems 
occur, has been significantly revised. Specifically, the field tables are 
now printed In a tabular format to provide easier location of spurious 
curvature or velocity values. Also, a specific print of the ERASE2 tenqjorary 
storage region is provided for each major calculation section of the program. 
All variables are identified by specific headings to locate critical para- 
meters. A full description of this output is given In Part II. 


2.3.3 Printer Plot 

Following the processing of the input boundary data and the construction 
of the Initial calculation grid, a one-page printer plot of the boundaries, 
primary orthogonals and streamlines is produced. Samples are given in Part 
II (Section 12.0). As indicated in these figures, the 5i~coordinates of the 
primary orthogonal lines and the of the dividing streamlines 

between channels are printed. The approximate locations of the boundaries 
are shown by plus signs. This allows convenient detection of gross irregu- 
larities in either the input boundary coordinates or the initial calculation 
grid. It also provides a map of the coordinate system which is a handy 
reference when examining the printout of the flow field solution. 


2.4 REVISIONS TO PROMOTE SOLUTION RELIABILITY 


To improve the starting reliability, three changes have been made in 
the computer code. First, the streamline positions are now cotqiuted at the 
zero-th refinement level. Second, the option for "constant density" calcu- 
lation has been enhanced so that it may be used when far-fleld and pressure 
boundaries are present. Finally, stagnation points are not treated as 
singularities (with adjacent ISTAG=3 points) until a well defined grid has 
been obtained in the stagnation region. 
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Previously, before any calculations were performed to correct the 
Initial guesses of the streamline positions, the grid was refined so that 
there would be at least one Internal streamline within each channel and one 
Intermediate calculation station between primary othogonals. (Primary 
orthogonals are defined In Section 3 of Reference 1). However, for cases 
where Inlet or splitter leading edges were closely spaced In the streanwise 
direction, it was found that the initial grid was not sufficiently smooth for 
the interpolation required. Therefore, the equations for the streamline 
positions are now satisfied before this first refinement. Irregularities 
in angle and curvature of the stagnation and trailing streamlines are now 
smoothed and a satisfactory base for the first grid refinement interpolation 
is provided. (For the case of a single channel, the initial streamlines will 
entirely coincide with the boundaries and the streamline corrections will be 
trivial.) 

Occasionally during the development of the grid in the early stages of 
refinement, an orthogonal will be located at a point on the boimdary where 
the curvature is very large. Also, unrealistically high curvatures will be 
found within the field because the zero-th grid is based on a flow/ area 
proportionality which causes kinks in the streamlines when an abrupt change 
in passage area (for exanqjle, a splitter vane) is encountered. When un- 
realistic curvatures exist over a sizable portion of the passage, the com- 
puted range in velocity will be very large and may give an indication of 
negative temperature (V>V TnaY ) or passage choking. To alleviate this possi- 
bility, the grid may be partially refined by assuming that the density every- 
where in the channel is equal to the stagnation density. With this assunq)- 
tion, choking cannot occur and (theoretically for incoiqjressible flow) V^ax - 
becomes infinity so a negative teii 5 )erature becomes inconsequential. This 
option is enforced by specifying a value of NODENS (see Part II) greater than 
(-1). For example, N0DENS=3 will call for the constant density assunq)tion 
to be utilized for the first three refinement levels. On the fourth refine- 
ment, the solution would revert to the full conq)ressible equations . At this 
point, the streamlines will be smooth and the resolution will be adequate 
to assure that the difficulties mentioned above will not be encountered. 

The default value of NODENS is zero. Therefore, the constant density 
approximation will be used for the zeroth grid solution unless the user 
specifies otherwise. This generally enhances the reliability of starting 
the STC solution except for high pressure ratio nozzle plume cases . In this 
case, the density level in the plume will be very different from the stagna- 
tion value and the use of the density option may retard convergence. There- 
fore, it is recommended that the default value of NODENS=0 be overridden 
by setting N0DENS=-1 when a nozzle plume is being calculated. 

The third feature now included in the program is to treat stagnation 
points as regular points until the grid becomes well defined. As such, a 
non-zero value of velocity will be conq)uted at "stagnation points" by using 
a numerically evaluated curvature defined by a three point fit. The point 
is converted to a true singularity point (where the streamline curvature is 
infinite and the velocity zero) when the numerically evaluated velocity 
falls to one half of the velocity at the adjacent point (on the same ortho- 
gonal). This adjacent point is then converted to an "ISTAG=3" point and its 
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position and curvature is then interpolated. The advantage of this new 
procedure is that it avoids the interpolation of such points when the grid is 
very coarse and, consequently, when the interpolation can lead to unreasonable 
values. Through the approximate solution of the streamline position equations, 
the new procedure "ties down" these previously interpolated points when the 
stagnation region is not well defined, thereby enhancing operational 
reliability . 


2.5 STRUT BLOCKAGE 


To account for the blockage effects of struts, additional terms have 
been included in the program logic. Specifically, for axisymmetric flow, the 
effective flow area is defined as follows: 

A = / 2irrXdS„ (14) 

C 


where S 2 is the distance along the orthogonal and X is the blockage factor. 
For planar flow, the above equation reduces to 


A = / XdS2 

C 


(15) 


Thus, in the axisymmetric case, the value of X is defined as the fraction of 
the unblocked circumference and is a value less than unity. In the planar 
case, X is the depth of the channel or flow stream in the third dimension. 
Values of X are input to the program defining the blockage at discrete points 
in a two-dimensional grid. The calculation grid point values are then 
obtained by linear interpolation. 


2.6 DESIGN OPTION-ARBITRARY PRESSURE BOUNDARIES 


In the version of the program described in Reference 1, provision was 
made for specifying a velocity along a far-field channel boundary. This 
capability allowed the numerically con^juted inner field to be matched to an 
economical, linearized far-field solution. In this case, the velocity from 
the far-field solution replaced the boundary condition that the last stream- 
line coincide with a fixed contour. In the present version of the program, 
these procedures have been generalized to allow user input of an arbitrarily 
specified pressure or velocity over a portion of any boundary contour. 

An approximate boundary shape must be specified to provide for the develop- 
ment of the initial grid. After several grid refinements, the program 
boundary conditions are switched to utilize the specified distribution of 
boundary pressures. 
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The pressure may be specified on up Co two separate boundaries . When a 
pressure boundary Is specified along one boundary, the flow balance procedure 
uses the specified pressure to set the velocity level. (This replaces the 
iteration for the velocity level generally required to meet the known channel 
cross-stream area). If boimdary pressures are specified along both channel 
boundaries the velocity on the upper boundary will dictate the cross-stream 
velocity levels (for the assumed curvatures) and an error will be computed in 
the lower boundary velocity. This velocity error is then Included on the 
right hand side of the correction equation (at the lower boundary point) and 
during the course of the iteration will be reduced to zero. In Appendix A, 
the streamline position correction equation for the lower boundary is formulated. 
The upper boundary correction equation Is the same as that used for far-fleld 
boundaries in the original version of the program. 

When using this option to design to a desired velocity distribution, it 
is advisable to first calculate the velocity distribution for an approximate 
(or initial) contour and then select a region over which the pressure will 
be modified. In the end portions of this region, the newly specified pres- 
sures should be faired smoothly into the previously calculated values. 

If this is not done, the pressures on the fixed portions of the contour will 
change abruptly. 


t 
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SECTION 3.0 


INTRODUCTION 


Part II of this report consists of the Users Manual for the STC Program. 
Detailed descriptions are given of the structure of the computer program, 
Streamtube Curvature Analysis (STC), and the program usage and operation. 

The information given herein represents an updated version of the program and 
supercedes that included in References 3 and 4. Two types of information are 
included in this User's Manual: (1) user oriented input sheets and output 

definitions along with an example case, and (2) programmer oriented descriptions 
showing program structure, program nomenclature, program messages and error 
codes, and operating instructions. The program listing appears as a separate 
document (Volume II) . 
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SECTION 4.0 


PROBLEM DESCRIPTION 


The Streamtube Curvature Analysis was formulated as a computer program 
to solve the invlscid equations of motion over a two-dimensional body (plane 
or axisymmetric) at transonic speeds. The computer program determines the 
flow field properties, streamlines, and pressure distribution over typical 
isolated nacelles and calculates the external pressure drag and the additive 
drag. (The additive drag is defined as the integral of pressure multiplied 
by the axial projection of the area taken along the entering streamtube from 
the undisturbed free stream conditions to the stagnation line on the cowl lip) . 
The solution is the direct type in that nacelle shape, mass flows, and flight 
Mach number shall- be the prime input data. A boundary layer procedure (SAB) 
is incorporated to allow evaluation of boudary layer displacement effects and 
friction drag as well as solution of the invlscid problems. 

The computer program is capable of analyzing the following geometries : 

a. Two-dimensional invlscid inlet problems without side spillage at 
zero or finite angle of attack. 

b. Axisymmetric zero-angle-of-attack isolated nacelle problem with: 

1. Short cowl nacelle in which fan duct air of a turbofan engine 
exhausts upstream of the primary air nozzle, and the pressure 
distribution on the aft nacelle (waist cowl) must be determined. 

2. Long duct nacelle in which the exhausts of both streams are 
confluent at the exit or mixed upstream at the exit. 

c. In all cases, the flow field may be calculated in the presence of 
a centerbody (or ranq)) whose leading edge may be positioned either 
forward or aft of the cowl lip plane. 

The flow field boundaries can be located as far upstream, and laterally 
displaced from the nacelle as far as practicable, to ensure minimum disturbances 
at the boundary and as far downstream as necessary to ensure correct nacelle 
trailing edge flow conditions. The program is capable of computing the inlet 
internal flow up to the assumed location of the engine face. The exhaust 
nozzle weight flow and aerothermodynamic properties of the exhaust flow shall 
be input quantities. The velocity of the exhaust may be sonic or greater. 

The computer program is structured to allow user control of the fineness 
of the computational mesh. This allows selective grid refinement in the 
stagnation line region and in regions where high velocity gradients occur. 

The most useful deck utilizes 768 grid intersection points and requires a 
central memory storage of 115,000 octal locations. 
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The size of central memory storage required is a function of the defined 
table sizes. These table sizes may be changed to meet the user's needs. The 
method to change the table size is included in this User's Manual. Otherwise, 
the basic logic of the STC program is identical and all the capabilities 
defined above are always available. 

The coding of the computer program meets the following requirements: 

1. The program has been written to run on the LRC's GDC 6400 and/or 
6600 computers, or any similar GDC 6400 and/or 6600. 

2. The bulk of the computer program has been coded in GDG FORTRAN 2.3 
language. Three subroutines have been coded in GDG Gompass 

1.1 language. The programs have been written to run under the SGOPE 

3.1 operating system. 

3. All input/output has been accomplished with GDG FORTRAN 2.3 state- 
ments. The standard system file names of INPUT for card reading 
and OUTPUT for printing have been used. In addition, input from 
tape files and output to tape files has been used. 

The description of the capabilities of STG do not include all possible features 
that are Included in the computer analysis. The user has control over the 
amount of grid refinement or computational mesh size, both by specifying 
local areas of mesh refinement and by setting the number of overall flow 
field refinements. The input geometry may be specified as coordinates only 
or coordinates and local surface angles on any boundary. Boundary layer 
effects may be selectively Included for any surface and imposed at any level 
of grid refinement. 
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SECTION 5.0 


METHOD OF SOLUTION 


5.1 BASIC EQUATIONS 

The following sections are concerned with the basic equations which 
are utilized in both, the STC "inviscid procedure" and the coupled turbulent 
boundary layer (SAB) procedure. 


5.1.1 STC 


The STC Program is designed to solve the equations of motion along stream- 
lines, T = constant lines, and along lines which are orthogonal to the stream- 
lines, C = constant lines. The variable C is introduced to avoid confusion 
with the velocity potential 4 > which is only applicable when the flow is 
irrotational. 


Across the streamlines, the continuity and Crocco form of the momentum 
equation are written: 


Continuity : 


Momentum : 


3A = ^ 

pv 


(C = Const.) (16) 


a = const.) (17) 


rjjj 9n 3n 


Along the streamlines the following forms of the energy and momentxam equations 
apply : 


Momentum: 
Energy : 


0 

DS 


0 

ds 


(T = Const.) 

(18) 

(1' = Const.) 

(19) 


where : 

A = Flow cross-sectional area = 2it r 9n 
C = Curvature of the streamline 

H = Stagnation enthaply 

n = Distance along the orthogonal 

p = Static pressure 

r = Radial coordinate 

S = Entropy 

s = Distance along the streamline 

T = Static temperature 

V = Velocity 
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4',W = Stream function, cumulative flow rate 

p = Density 

The solution method is an extension of the conventional streamline curva- 
ture method. It may be briefly described as follows: First a crude grid of 

streamlines and orthogonal lines are assumed. (Please refer to Fig. 19 , 

Second, the curvature of the streamlines at each of the grid points is evaluated 
Third, the momentum equation is Integrated along a line normal to the stream- 
lines to obtain velocity and the continuity equation is integrated to determine 
the "correct" streamline positions (for the assumed curvature field) . These 
are indicated by the "x" in Fig. 19. Fourth, an adjustment, 9n, is computed 
by considering ( 1 ) the difference between the computed and assxmed streamline 
positions and ( 2 ) the effect of the implied curvature modification in the 
Integrated momentum equation. Finally, the streamlines are repositioned by 
the 9n values. 

Because the movement of any one grid point alters, through a change in 
curvature, the velocity at nearby points, it is hl^ly desirable to account 
for these interrelating point adjustments simultaneously. The utilization 
of a simultaneous solution procedure, employed here, is not part of the 
classical streamline curvature method [5,6,7]. In comparison, the classical 
method yields calculation times which are very slow, especially for a closely 
spaced calculation grid. In concept, the set of simultaneous equations for 
the normal streamline adjustments are formulated from the finite difference 
equivalent of the following equations: 


a^On) . (1-M^) a^On) 

2^2 2 

. (pv) as 


( 20 ) 


where : 

8n = 

s = 
M 

pV = 
F 


Required streamline adjustment in the normal direction 
Stream function 

Curvilinear distance along a given streamline 

Mach number 

Flow per xmlt area 

Driving (or error) function derived from the solution to the 
Integral continuity and normal momentum equations. 


This equation is derived in (Ref. 1) for the special case of isentroplc 
2-ditt»nsional flow. (These limiting assrnnptions are utilized only to main- 
tain simplicity of illustration; they are not part of the computer program) . 
From a mathematical point of view, the above equation is similar to the con- 
ventional equations for velocity potential or stream function, namely. 
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2. Evaluate Curvature 


3. Integrate The Cross- 
stream Momentum Eqn. 
and Continuity Eqn. To 
Determine Hie "Correct" 
Streamline Positions. 



4. Solve The Matrix 

Equation for 6n and 
Move Hie Grid Points. 


Figure 19. Solution Technique. 
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( 21 ) 




* y 


dy 


d\ 

3x^ 


= 0 


( 22 ) 


X and y are coordinates in a two-dimensional rectangular system and <t) and '1' 
are the velocity potential and stream function, respectively. (Again, all of 
the equations have been restricted to Isentropic flow for illustration) . 

However, for the purpose of calculating transonic flow, the use of Eq. 

(5) offers a distinct advantage over Eqs. (21) or (22) for the following reason: 
Because the grid is always aligned in the streamline and normal to streamline 
directions, no cross-derivatives appear. In consequence, the finite difference 
star is simply switched from a subsonic representation to a supersonic 
representation Illustrated in Fig. 20. 

Notice that for supersonic flow, no points downstream of the orthogonal 
line are included. This reflects the physical reality that disturbances dovm- 
st re am will not affect the flow upstream. It is, of course, because the 
coefficient term, (1-M^) , passes through zero and changes sign that the star- 
switching noted in Fig. 20 is appropriate. 

If the grid system is not aligned with the flow direction, a cross 
derivative appears in the equation, as in Eq.(20). Unfortunately, the mixed 
derivative coefficient, does not have the same sign change property 

and, therefore, star-switching with this equation is more difficult for 
general flows with some angularity. 

The extended streamline curvature method, here referred to as the Stream- 
tube Curvature (STC) method, has then the following features. 

• No additional complexities arise when the flow is rotational. 

• The slip lines between the exhaust jet and the external flows can 
be handled precisely. (The procedure is to consider two coincident 
streamlines. Their position and pressure are the same; their ve- 
locity and stagnation properties may be different) . 

• From a numerical point of view, the streamline/orthogonal line 
oriented grid facilitates the analysis of transonic fields, as 
described above. 

• The streamline/orthogonal line grid also provides a mapping of the 
flow field into a rectangular domain. This is helpful from the 
standpoint of coin)uter program organization. 
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— ► M <1 


— ^ M = 1 



Figure 20. Finite Difference Stars for Subsonic 
and Supersonic Flow. 
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• Boundary layer displacement effects are easily included by coupling 
the STC Inviscid solution with a turbulent boundary layer procedure. 

The STC Program has also been designed 

• to handle multiple streams , 

• to adjust the flow rates of the jet exhaust streams to meet the 
"Kutta" condition at trailing edges or the (2-dimensional) choking 
condition, 

• to place grid points at locations in the flow field where they are 
needed, as determined by local variations of the dependent variables, 
and 

• to allow external flow analysis by incorporating matched near-field 
and far-fleld solutions. The far-field solutions are obtained 
analytically, utilizing small perturbation theory. 

‘’•1-2 SAB - (Turbulent Boundary Layer ) 

The boundary layer procedure cotq)led with the STC inviscid solution is 
the method of Stratford and Beavers (SAB) described in detail in References 1 
and 2, In the SAB method, the Integral bomdary layer parameters 6, 6*, 6 
are expressed in closed form as a fvinctlon of Mach number, equivalent flat 
plate length, and Reynolds number based on the equivalent flat plate length. 


h 

6* = f(M) X Re" 

.<5 . 


C23) 


The equivalent flat plate length in the above expression is defined as 
the length over which a boundary layer growing on a flat plate at the given 
Mach number would acquire the same thickness as the real boundary layer 
at that given location. For axisymmetric or plane flow 



where = Sw = 


distance measured 
along the wall 


a 

a 


Axisymmetric 
0 Plane 


(24) 


P 
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The specific working formulas for calculation of the integral parameters 
are as follows; 


For 


10^$Ry$10^ 

0 = 0.036 (I+.IM^) 


-0.7 -1/5 

X 


0.44 -1/5 

6* = 0.046 (1+.81T) X R^. 

6 = 0.37 X 


For 10^^Rx<10^ 


-0.7 -1/6 

0 = 0.022 (l+.llT) X R^ 




(25) 

(26) 

(27) 

(28) 


6* = 0.028 


(1+.8M^) 


0.44 -1/6 


X R- 


X 


(29) 


6 = 0.23 X Rx"^'^^ (30) 

With these relations, the distributions of 0(S ) and 6*(S ) can be 
calculated, given the boundary layer edge pressure or velocity'’'^distribution 
form STC. The local skin friction coefficient for determination of the 
friction drag is evaluated numerically using the integral momentum equation. 


C. = 2 


d0 

dS 


w 


I 


dV 

dS 


( 26 H- 6 *) + 


w 


si 

f 


dr 

dS 


w 


2^ ^ 

p dS 


w 


(31) 


where e 


0 Plane flow 

1 Axisymmetric flow 


Boundary layer separation is detected during the calculation by evalua- 
tion of the Stratford separation parameter (Reference 23) . This quantity is | 
defined as : 



V dS 


w -J 


1/2 



where Cp = 1 



(32) 


and is evaluated at the minimum pressure 
point of an adverse pressure gradient region 
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The Stratford parameter is only calculated in a region of adverse pressure 
gradient (dP /dS )>0. The disttince along the surface is taken as the distance 
from the begYnniSg of the adverse pressure gradient region, and M, is evaluated 
at the minimum pressure point. For practical purposes, separation is assumed 
to occur if F attains a value of 0.5 or greater. The calculation continues 
beyond the separation point, but the separation point and all values downstream 
are flagged in the boundary layer output and the STC output if the problem is 
restarted to include displacement effects. 


5.2 AN OUTLINE OF THE CALCULATION SETUPS 

The operations performed by the STC Program may be outlined as follows: 

1. Define the flow regions and locate (approximately) the "primary" 
orthogonals and the streamlines which divide the internal and 
external flows. 

2. Refine the grid as required by Inserting additional streamlines and 
orthogonal lines between those already existing. 

3. Compute the streamline angles and curvatures. 

4. Compute the orthogonal line angles and move the grid points along 
the streamlines to obtain orthogonality. 

5. Compute the velocities on the "far-field" bomdary. 

6. Adjust the flow rates in the e::diaust streams, if any, to meet the 
calculated choking flow rate. 

7. Integrate along each orthogonal the momentum and continuity 
equations [Eqs. (16) and (l7)]. 

8. If the streamline positions are within a "rough tolerance", adjust 
the flow rate of variable flow channels to meet the "Kutta" condi- 
tion at trailing edges. Step 7 is repeated with each iterative 
adjustment of the flow rates. 

9. If the streamlines are within tolerance, return to Step 2 for 
additional grid refinement (unless grid refinement limits have 
already been reached). Otherwise, continue to Step 10. 

10. Determine if the streamline positions are within final tolerance. 

If so, jump to Step 14. Otherwise continue to Step 11. 

11. Set up the matrix equation for the streamline correction, 6n. 

12. Solve the matrix equation. 
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13. 


Modify the streamline position by 6n, and return to Step 3. 

14. Calculate and print the output quantities; calculate boimdary 

layers and adjust. wake table at trailing edges; then return to 
Step 1 for the next case, if any. 

The first operation includes reading the card input for a description of 
the geometry and flow properties. The computer program has been written to 
have general capability for analyzing a great variety of configurations . The 
first step in the programmed logic is to develop a table of orthogonals or 
calculation stations for each of the several flow regions. The regions are 
determined as illustrated in Fig. 21 so the calculation can proceed from 
upstream to doxmstream. The boundary of each region is defined as a primary 
orthogonal. As shown in Fig. 22, the initial grid which is developed con- 
tains only the primary orthogonals and the double streamlines which separate 
the various streams. A full inner Iteration is carried out at this point 
to improve solution reliability at later stages of grid refinement. 

The second step in the computational procedure is the grid refinement. 

The very crude grid, obtained in Step 1, is refined before the second solution 
of the flow field equations is executted. A new orthogonal is placed within 
each region and, likewise, a streamline is inserted in the middle of each 
channel. In the external channel, additional streamlines are placed close 
to the body. After the solution has been obtained for this net, the grid 
Intervals are halved as required. This may be likened to the steps taken 
when one "flux plots" a flow field by hand. First, major flow lines and 
normals are sketched in, and then more and more streamlines and orthogonal 
lines are added until the desired resolution is obtained. At each step in the 
process the positions of the lines are adjusted to meet the correct solution 
requirements. The procedure automatically provides for grid refinement in 
regions of higb curvature and high acceleration or deceleration. The stream- 
line and orthogonal lines which are added between existing lines are not 
required to span the field if only local refinement, near the body, is required 
The refinement procedure presently built into the program uses a criteria 
involving the distance and velocity increment between grid points. These 
refinement criteria are discussed in detail in Sections 11 and 13. 

The third step in the method is to determine the angles and curvatures 
of streamlines at each grid point. For subsonic portions of the flow field 
this is performed by fitting a piecewise continuous cubic polynominal in a 
coordinate system which is locally rotated for each interval. The resulting 
fit is analogous to the curve produced by a beam which is loaded by discrete 
forces so as to pass through the given grid points. The locally rotated co- 
ordinate system removes the restriction that requires the slope to be small. 

For grid points located in a supersonic region, backward difference formulas 
are employed. Either 3-point, 4-polnt, or 5-point formulas may be optionally 
selected. Again the coordinate system is rotated so that slopes in the curve 
fitting coordinate system are small. 
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Figure 21. Subdivision of the Flow Field into Regions, 
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Figure 22. Initial Streamline/Orthogonal Grid Before the First Refinement 



In the fourth step the orthogonality of the grid points Is checked and 
points are moved along the spline curve as required to achieve normal Inter- 
sections between the two sets of lines. Also, the normal distance, n, is 
computed for each grid point as measured from the lower boimdary of the 
orthogonal. If boundary layer data are available from a preceedlng STC pass, 
displacement effects are included during the orthogonalization procedure. 

Interpolated corrections are applied to the boundaries using the following 
^ relations ; 



Rc = Rb “ 

Zc = Zjj + 6* sin (<!>(,) 


) 

Lower Boundary 



Rc = Rb + 5* cos (<{ic) 
Z(, = Zj, - 6* sin (<|)c) 


If the boundary layer on the solid surface is separated, a comment to this 
effect is printed each time the boundary is accessed to calculate an orthogonal 
intersection. When this situation is encotmtered, the user is advised to dis- 
continue the calculation, since the displacement thickness Information down- 
stream of the separation point is in error. 
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When the initial grid is set up, a boundary is placed some distance 
avay from the body. This boundary becomes the interface between the near- 
field and far-field solutions. The near-field is computed by the stream- 
tube curvature method and the far-field is cor 5 >uted by linear small perturba- 
tion theory. In the process of Iterating, this boundary (which is also a 
streamline) will float so that its shape and velocity distribution are matched 
by both the inner and outer solutions. In practice, the shape of the inter- 
face streamline (also referred to as the far-field botmdary) is first assumed. 
Using the far-field equations, the velocity distribution is calculated. This 
is Step 5. These velocities are subsequently employed in the near-field 
analysis and from this comes a revised shape for the interface streamline. 
Revised velocities will then be conqmted in Step 5 during the following 
iteration cycle, and so forth. 

Step 6 is the modification, as required, of the flow rates of the exhaust 
streams. For boattall analysis of nacelles, the internal geometry of the 
exhaust passage is required input to the STC program. Because of streamline 
curvature effects, the discharge coefficient for the nozzle will be somewhat 
less than unity. The user, however, may input a flow rate based on unity 
discharge flow coefficient or, for the matter, any approximate value. De- 
termination of the velocity distribution across the throat of the nozzle 
will be determined within the STC framework and the evaluation of the maximum 
"choked" flow rate is Step 6 of the calculation procedure. 

Step 7 is the solution of the flow field equations per se. This section 
of the program is referred to as the "flow balance"; Eqs . (.16) and (17) are 
Integrated. In the external regions of the field, the momentim equation is 
integrated from the far-field interface boundary to the body (or to the 
centerline or lower boundary, whichever exists). The integral form of the 
momentum equation is : 


An vj- 


An Vpp 


2 



(33) 


where : 

Vpp = Velocity as determined in Step 5 along the far-field streamline. 

Vq = Velocity at any streamline with orthogonal distance nQ. 

n^j, = Distance measured along the orthogonal to the far-field 
streamline. 
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Although not reflected in Eq. (33), the effect of varying total pressure 
behind a shock wave is also Included. Also, if a slip-line occurs in the 
field, the velocity jump equations. 



(34b) 


(34c) 

are employed where the subscripts (+) and (-) denote conditions on the stream- 
lines above and below the slip-line, respectively. 

The velocity, total temperature and total pressure allow determination 
of the density at each grid point, and the Inverse product of density and 
velocity is integrated to find flow area. 


V 


V 2 

Cp =H_ 



(35) 


Step 8 consists of the iterative flow adjustment of the variable flow 
channels to meet the departure streamline pressure compatibility condition 
("Kutta" condition) at all trailing edges. Step 8 is executed only if the 
streamlines are within tolerance. Also, Step 7 is repeated with each 
iterative ad j iis tment of the flow rate. 

The cumulative flow areas calculated by Eq. 35 are compared with the 
geometric areas of the streamlines used in Step 3. The difference between 
these two values is used as a convergence check (Steps 9 and 10) and in the 
streamline correction equation. Step 11. 

In Steps 11 and 12 the proper adjustment of the streamline positions is 
determined and in Step 13 the grid points are moved in the normal direction 
by this computed adjustment. 

The iterative sequence is to start with a crude grid, as noted above, 

^d to repeat Step 3 and througji 13 until the flow balance error is small. 
This is often accomplished in one or two iterations . Any variable weight 
flows are then converged. The grid is then refined to the next level and the 
field is reconverged. The refinement /convergence process is continued until 
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the grid refinement criteria are satisfied, or alternately, until computer 
storage limits are reached. At this point, additional loops through Steps 3 
to 13 may be performed until the flow balance error is satisfactory. 

In Step 14, output quantities are calculated and printed. Boundary 
layers are also calculated and stored for a subsequent STC restart. Finally,, 
the wake table is adjusted to reflect boundary layer displacement effects at 
trailing edges . 

A complete description of the details of the numerical procedure may be 
found in References 1 and 3 through 9. 
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SECTION 6.0 
PROGRAM DESCRIPTION 

6.1 DATA STORAGE TABLES 

The framework of the STC program is designed to allow flexibility as to 
the configuration to be analyzed. For example, very weak limits are placed 
on the number of flow boundaries and the mnnber of channels into which the 
flow is split. And no specific limits are placed on the ntimber of stream- 
lines or the number of orthogonal lines in any given region of the flow. 

To accomplish this, the bulk of the calculation data is saved in arrays 
which are singly dimensioned. Within each array the data are packed together 
to maximize storage efficiency. Descriptions of these tables are provided 
below. 

6.1.1 Field Tables 

Flow field data are stored in singly subscripted arrays, one for each 
variable. The quantities saved at each streamline - orthogonal line grid 
intersection are as follows : 


Z 

R 

VM 

PHIl 

CURV 

51 

52 
B 

RHS 

DS2 


axial position, abscissa 

radial or vertical position, ordinate 

velocity 

streamline angle, measured from horizontal 
curvature of the streamline, -d(l>/dS]^ 
curvilinear distance along the streamline 
curvilinear distance along the orthogonal line 
coefficient in matrix equation for DS2, Indicator 
of subsonic (B>0) or supersonic (B<0) velocity 
right hand side of matrix equation for DS2 
correction of streamline positions 


In the Fortran coding M is the symbol commonly xised as the subscript for the 
above arrays and NM is the number of field points (or maximum value of M) . 
The points are grouped by orthogonals starting at the upstream orthogonal. 
Points along the orthogonal are grotjped together. Hence, the point below 
and above the M^b point are (M-1) and (M+1) respectively. The neighbors in 
the streanwlse direction, however, are determined by referring to the JMS- 
table described in Section 6>2. 
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6.1.2 Channel Input Data Table 


Information such as the boimdary coordinates and flow properties is 
compactly stored in a single array, TABLES, so that only the total amount of 
information saved is. limited by the array size. No limit is placed on the 
amount of information to be placed in any one Table of which there are seven; 

• channel input table 

• boundary table 

• table of convected properties 

• table of wake displacment thickness 

• boundary layer data table 

• flow adjustment table 

• station table 

In the first table, CHDATA, input information read from page STC/sheet-3 
is stored. The information is stored in subtables, one subtable for each 
channel, and the arrangement of each subtable is as follows: 

CHNAM channel name (BCD) 

LHNEXT length of subtable 

WTFLOW flow rate (if input) 

TTO total temperature 

PTO total pressure 

TSO static temperature 

PSO static pressure 

MACHO Mach number 

AO flow cross-section area 

VARY indicator as to whether the flow rate may be varied 

RG gas constant 

GAM ratio of specific heats 

NR not operational 

NC not operational 

TAB not operational 

BB(NR&NC) not operational, zero length array 
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Except for the first two words, CHNAM and LHNEXT, all of the above Input 
items are optional and they are all equal to BITS unless values are supplied 
according to input instructions for STC/Sheet-3. If data for a second channel 
are supplied, these data will follow the first channel and so forth for any 
number of channels. The first word of the CHDATA table is at location LHO 
and the last word is at location LHE relative to the origin of the TABLES- 
array. If LHE=LH0-1, no channel data has been input and the channel data 
table has zero length. In the Fortran coding, the subscript LH is used to 
refer to the channel data. 

Channel information is read into the STC program by subroutine RCD; the 
stored channel data is utilized by routines RTCFI, BCONV, ADJWF, ISBOT, and 
SABBL. 


6.1.3 Boimdary Table 


Directly following the channel data are the coordinates of the boundaries 
Again a subtable for each boundary is constructed and the information is 
stored in the following order: 


BDT 

LBNEXT 

LBZl 

CHNAME 

UP 


LEDEX 

ZBT(l) 

RBT(l) 

ANGBT(l) 

ZBT(2) 

RBT(2) 

ANGBT(2) 

o 

o 

o 

o 


boundary name 
length of subtable 

index increment to first coordinate in ZBT, RBT, ANGBT-lists 

channel with which the boundary data is associated 

upper or lower boundary indicator (if the boundaries aroimd 

a leading edge have been collated together, then CHNAME(2)=UP 

is the name of the channel above the leading edge and 

CHNAME (1) is the name of the channel below the leading edge) . 

index (relative to ZBT, RBT, ANGBT) of the leading edge point 

when boundaries are collated 

axial coordinate (x) 

axial coordinate (y) 

surface angle (measured from x-axis in radians) 
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BDNAME 


name of a specific boundary when several boundaries are 
collated together in one subtable 
LBA, LBB Index limits relative to ZBT, RBT, ANGBT of the 

coordinates for BDNAME 

(NOTE - the above 3 items are repeated for each boundary where LBZl/3 
is the number of collated boundaries. The existence of this information 
results in the displacement of the ZBT, RBT, ANGBT coordinates to hi^er 
memory locations. BDNAME, LBA and LBB are equlvalenced to ZBT, RBT 
and ANGBT respectively) . 

Boundary coordinates are necessary input to the program so two or more 
subtables will always exist, stacked one after another. The first word of 
the BDYTAB table is at location LBDO and the last word is at location LBDE. 

LB is the subscript vised to refer to the boundary data. 

Boimdaiy coordinates as supplied on page STC/Sheet-2 of the input sheets 
are read by routine RBD. This routine converts the angles from degrees to 
radians, translates and rotates the coordinates as required and stores the 
points so that as one proceeds from one point to the next (or walks along the 
contour) one's left arm is next to the flow field. Thus, the ordering of 
points on an upper flow boundary will be reversed and the angles incremented 
by 180 degrees. 

The several boundaries which may be defined to comprise one continuous 
contour are collated in subroutine BLDTAB. The boundary table is referenced 
in routines BDYPTM, LBDYBL, and BLTBBL. 

6.1.4 Table of Convected Properties 

For each channel, a subtable of convected properties and channel flow 
data is built by subroutine BCONV from data in the channel input table (if it 
exists) and input data from STC/sheet 1. This table contains some of the 
same information as the CHDATA table. Information contained in the CONVTB 
table is complete and follows a consistent format. 


CH 

channel name 

LTNEXT 

length of subtable 

NPT 

not operational (=1) 

LPSI 

(=15) 

LTT 

(=16) 

LPT 

(=17) 

LRCU 

(=18) 

CRG 

gas constant 

CPGJ 

specific heat at constant pressure 
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C2CP 

twice CPGJ 

QGAM 

inverse of ratio of specific heats 

FGT 

(y-i)/y 

FGP 

y/(y-i) 

FGR 

i/(y-i) 

AREATB 

area for calculating flow rate 

PSI 

flow rate for the channel 

TT 

stagnation temperature 

PT 

stagnation pressure 

RCU 

not operational 


Again the data ,for the several channels are stacked one after another. LTO 
and LTE are the first and last locations of the table relative to the origin 
of TABLES. LT is the subscript used to retrieve convected property informa- 
tion of sxib routines BLDTBS, MCONV, ADJWF, TTPT, and SABBL. 


6.1.5 Table of Wake Displacement Thickness 


The wake displacement table is constructed if there exists a trailing 
edge , and is arranged as follows : 


X2W 

LWNEXT 

SIW 

(1) 

SIW 

(2) 

SIW 

(N) 

DST 

(1) 

DST 

(2) 

DST 

(N) 


streamline coordinate, C 2 

length of subtable, = 2N+2 

list of distances from trailing edge 


list of wake displacement thicknesses 


This table is built in the subroutine BLDTBS by a call to subroutine 
BWAKE. DST is determined so that the wake thickness equals the trailing 
edge thickness at S1W=0. Thereafter the wake thickness decrease at the rate 
of 0.1 times the distance from the trailing edge. At 10 trailing edge 
thicknesses downstream, the wake thickness is zero. 


As many wake thickness subtables are built as there are trailing edges . 
The wake thickness table begins within the TABLES-array at location LWO 
and ends at LWE. The table is referenced by the TTPT routine. 
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At the completion of a bovmdary layer calculation (if specified) , the 
wake table entries at the trailing edges are adjusted to reflect the trailing 
edge displacement effects. This procedure is accomplished by subroutine 
RBWAKE. 


6.1.6 Botmdary Layer Data Tables 

Bovmdary layer data in the STC-SAB program are stored in the TABLES 
region immediately before the flow adjustment table. A sub table for each 
boundary layer is constructed and the Information is stored in the following 
order: 


BNAME 

LBLNXT 

NSEP 

DUMMY 

SWEEP 

SIGN 


SW(1) 

DSTAR(l) 

DDSTAR(l) 

SW(2) 

DSTAR(2) 

DDSTAR(2) 


Boundary name 

Pointer to the next boundary layer table 

Index pointing to separation location (normally 0) 

Reference distance for alteration of coordinates in the 

bovmdary layer table; viz, boundary origin 

Boundary type 

-1 Upper boundary 

+1 Lower boundary 

Distance along surface 

Smoothed displacement thickness (6*) 

Slope of smoothed displacement thickness (d6*/ds) 


The boundary layer data table is located between LDO and LDE in the 
tables region. The index limits LDO and LDO are stored in /IXORIG/af ter 
LSE. LD is the subscript used to reference the bovmdary layer data tables 
(subroutines BLTBBL and BDYPTM) . 


6.1.7 Flow Adjustment Table 

The flow adjustment table is also created in the BLDTBS routine, one 
subtable for each trailing edge. The information contained is as follows: 

XIF orthogonal coordinate of the t.e., Cl 

X2F streamline coordinate of the t.e., ^2 

XIBF Cl coordinate of the choked station of the flow above the 
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t.e. if not XIF 



XIAF 


SIF 


NCHB 

NCHA 

JORDER 


VNR (12) 


Cl coordinate of the choked station of the flow above the 
t.e. if not XIF 

curvilinear streamline distance to the trailing edge (along 
the upper surface of the body) . This value is used for 
interpolating the wake displacement thickness, 
number of channels below the t.e. 
number of channels above the t.e. 

= -1 if the single channel flow is choked 

= 1 if channel flow rates below the t.e. are known 

= 2 if channel flow rates above the t.e. are known 

12 element storage array used by subroutine NEWRAP for the 
flow iteration 


The data stored in this table is used by subroutine ADJWF which adjusts the 
flow so that at each trailing edge the pressure difference from one side 
of the trailing edge to the other is reduced to zero (to satisfy the Kutta 
condition) . If this condition cannot be satisfied, the flow in one of the 
channels will be choked and routine ADJWF adjusts the flow to the maximum 
choked value. The flow adjustment table is located between LFO and LFE. 
Each subtable is NFCOLS (=20) in length. 

6.1.8 Station Table 


The station table is the last of the compacted tables and it contains 
information for each orthogonal line. The table grows during the calculation 
process because the number of orthogonals is increased to obtain a refined 
grid. Because it is the last table it can easily be extended into the unused 
portion of the allotted memory. The data saved for each of the orthogonals 
is arranged in the station table as follows: 


XI 

LNEXT 

MLB 

MQB 

PRIM 


station coordinate , 5^ 
length of subtable 

field point index of first (lower boundary) point of the 
orthogonal 

field point index of the last (upper boxindary) point of the 
orthogonal 

primary station indicator, T or F, (a primary station, is one 
of the original grid stations and it will pass through boundary 
end points) . 
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TYPELB 

NAMELB 

ILB 

FLB 

SILB 

TYPEUB 

NAMEUB 

lUB 

FUB 

SlUB 

VMB 

DWDV,S CHOKE 


X2CL 

SLSWI 

MCL 


type of the lower boundary, l.e. SOLID, FIELD, FARFLD, etc. 
for indicating proper boundary condition 

name of the lower boundary used for referencing the boundary 
table 

boxmdary (table) Interval of the orthogonal boundary 
intersection 

fractional position in that interval 

curvilinear distance from beginning of the interval 

type of the tipper boundary 

name of the upper boimdary 

Interval of the upper boundary 

fractional position in the Interval 

curvilinear distance from the beginning of the interval 
boundary velocity, used as an initial guess for the velocity 
iteration in the FLOBAL routine 

slope of the flow rate versus velocity curve used for the 
velocity iteration in the FLOBAL routine, = SCHOKE if flow is 
choked at this station 

^2 coordinate of the control streamline used for positioning 
orthogonal lines in the PTMOVE routine 

sonic line-shock wave indicator, = 1.0 for mixed flow, = 0. 
otherwise 

field point index of control streamline, used in PTMOVE only 


This group of 20 items is repeated for each station, starting with the 
upstream and proceeding to the downstream stations. Station table entries 
at trailing edges contain seven additional items; 


ANGTE Boundary surface angle 

PTTE Total pressure 

PSTE Static pressure 

FGRTE Function of y, 1/(1-Y) 

RGTE Gas constant 

ANGEXP Flow departure angle just downstream of the trailing edge 
BSQEXP M^-1 just downstream of trailing edge 
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Quantities above the trailing edge are stored with the station above the 
trailing edge and quantities below the trailing edge are stored with the 
station below the trailing edge. At leading edges, the station table con- 
tains 22 entries. The additional items in this case are; 

CITRVLE Leading edge curvature 

ANGLE Leading edge angle 

Field point calculations are performed by looping through the station 
table, starting with the first station at L=LO and proceeding to the last 
station. The last word of the last station is at location LESTA. With excep- 
tion of the boundary layer table, the above tables are initially constructed 
with gaps between them. The channel data and boundary tables are constructed 
simultaneously in the ROD and RDC routines called from REDINP. The ma xi mi im 
size of the channel data region is preset by the value of MAXLH. MAXLH is 
Initialized to 400 but may be input as any value if such should be necessary. 

After reading the card data input, at the beginning of the BLDTAB routine, 

the boundary data is moved down in memory so that it is just above the channel 
data. The gap between the two tables is thus eliminated. 

The remaining four tables are then all built in subroutine BLDTBS. 

Again spaces between the tables are provided so that they can be constructed 
simultaneously and grow to the following lengths without interfering with 
each other: 



Input Variable 

Default 

Table Name 

Name 

Value 

Convected Properties 

MAXLT 

200 

Wake Displacement Thickness, 

MAXLW 

200 

Flow Adjustment 

MAXLF 

200 


At the end of the BLDTBS routine the gaps between the tables are 
eliminated by moving each table down in memory. The result is that all 
information contained in these tables is stored compactly, and the last 
table, the station table, has space into which it may expand. The length 
of all seven tables combined is limited to the length of the TABLES-array, 
and this value is adjusted at program load time to meet the size requirements 
of the flow field. 


6.1.9 Streamline Table 

The streamline table is not stored in a compact arrangement. It consists 
of three arrays in a labled common SLTAB. 

W cumulative flow rate for the channel 

X2 ^2 ~ coordinate 

SLCHN channel name of which this streamline is a part 



The streamline number, J, is the subscript used to access information 
in each of these three arrays. All of the streamlines for a given channel 
are together and in order (proceeding away from the centerline) . However, 
no special ordering of channels is required. 


6.1.10 Table of Leading Edge and Trailing Edge Points 

The Leading Edge-Trailing Edge Point (LETEPT) table is constructed in 
routine BLDTAB and used in routine BLDTBS for the purpose of defining the 
flow regions and primary orthogonals. The data in this table is not saved 
after the BLDTBS routine is executed. 

The information in this table is obtained from the boundary table. The 
ends of boundaries and double points contained within the boundaries are 
listed, together with the boundary and channel names. Each "line" of the 
LETEPT table contains the following ten items. 

axial coordinate of boundary end point of double point 
radial or vertical coordlnr.ee 
mean angle of the flow at XE, YE 

number of times the same point has occurred as an upstream 
boundary end point. Normally NLE=0 or 1. If NLE=2 then the 
point is a leading edge. 

number of times point has occurred as a downstream end point. 
NTE=2 for trailing edge point. NLE=NTE=0 for a double point 
in the boundary table 
name of flow stream above point 
name of flow stream below point 
name of bovmdary (UPPER=F) 
name of boundary (UPPER=T) 

number of times the point has been used in developing the 
ORTCHN-table, initially = -1 for a double point 

After the table is constructed the points (or lines) are sorted so that 
the upstream points are first and the points follow from upstream to down- 
stream in order. If difficulty is encountered in the development of the 
basic grid, the information contained in the LETEPT-table may be helpful 
in diagnosing the error. 


XE 

YE 

ANGE 

NLE 


NTE 


CHL 

CHU 

BDL 

BDU 

NUSED 
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6.1.11 Table of Channels Embraced by Each Orthogonal 


The LETEPT-table contains all points through which (primary) orthogonal 
lines are to pass. From this table the ORTCHN-table Is developed as an aid 
to the construction of the Initial grid of streamlines and orthogonal lines . 
The latter table contains a list of all of the channels embraced by each 
orthogonal. Specifically each "line" of the table contains: 


LEDGE 

index of the point in 

the LETEPT-table 

LRPREV 

previous line number 
orthogonal 

(or ORTCHN-table index) or the upstream 

CHNA(l) 

channel names 


CHNA(2) 

chsuinel nait^s 


CHNA(3) 

channel names 


CHNA(A) 

channel names 



The above Information Is tabulated (one line) for each orthogonal. Also, 
the first two lines are developed as dimmies for the purpose of listing all 
channels. Other dummy lines may also be Inserted. LRD Is the Index Incre- 
ment between lines and Is equal to the total number of channels minus two. 


6.1.12 Table of Index Limits 


The labeled common /IX0RIG/ contains the Index limits for each of the 
above listed tables. The Items and order of storage are as follows: 


LHO 

location 

of 

LHE 

location 

of 

LBDO 

location 

of 

LBDE 

location 

of 

LTO 

location 

of 

LTE 

location 

of 

LWO 

location 

of 

LWE 

location 

of 

LFO 

location 

of 

LFE 

location 

of 

LO 

location 

of 

LESTA 

location 

of 

LSO 

location 

of 

LSE 

location 

of 


channel table origin 

channel table end 

boundary table origin 

boundary table end 

convected properties table origin 

convected properties table end 

wake displacement thickness table origin 

wake displacement thickness table end 

flow adjustment table origin 

flow adjustment table end 

the station table origin 

the station table end 

the shock point table origin (not used) 
the shock point table end (not used) 
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LDO 


location of boundary layer table origin 

LDE location of boundary layer table end 

LDUM(4) unused 

MO unused 

I 

NM number of field points 

NJ number of streamlines 

NFCOLS number of "columns" in the flow adjustment table 
MAXNJ maximum number of strealines (dimensional limits are 
MAXNJ=128) 

MAXOL maximum number of points on any one orthogonal (dimensional 
limits are MAX0L=96) 

MAXNM maximum length of field arrays, calculated in subroutine 
REDINP 

MAXLE maximum table length, calculated in subroutine REDINP 

LEO location of the first word in the LETEPT-table , (=1) 

LEE .location of the last word in the LETEPT-table 

\ 

LRO location of the first word in the ORTCHN-table 

LRE location of the last word in the ORTCHN-table 

LRD the ORTCHN-table is subdivided into lines and LRD is the length 

of the lines 

6.1.13 Boundary Layer Input Table 

The boundary layer input table is stored in labeled common BLBDY*. Input 
boundary layer information as supplied on page STC/Sheet-2 of the input 
sheets are normally read by routine RBD. The resulting table consists of 
the following three items stored serially for each boundary: 

BLB(I) Boundary name. 

BLB(I + 1) Indicator designator whether a boundary layer calculation 

is to be performed. 

0 - No 

1 - Yes 

BLB(I + 2) Initial equlvzilent flat plate distance to first point 
on boundary. 
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The subscript I ranges from 1-58 and is incremented by 3 for each boundary. 
Information for a maximum of twenty (20) boundary layers may be stored; viz, 

Common/BDBDY/BLB(60) 

t 

The data input for this table may also be read by subroutine REDINP if 
the boundary layer is to be initiated by restarting a "non-boundary layer" 
run. 


6.2 ACCESSING DATA IN THE FIELD TABLES 

The data stored in the field tables consist of the Z, R, SI, S2, PHIl, 
CURV, VM, B, RHS, and DS2 at each point in the field. An additional array 
termed the JMS table provides access to the information in the field tables. 
This information (subscripted M=1,NM) consists of: 

Streamline number 

M subscript of upstream grid point; MU = 0 at beginning of 
streamline. 

M subscript of downstream grid point; MD = 0 at end of 
streamline. 

Point type Indicator 
ISTAG = 0 - Normal point. 

1 - Stagnation (or singularity) point. 

2 - Trailing edge point or a point fixed on the 
body surface used to locate a primary 
orthogonal. 

3 - Point adjacent to a stagnation (singularity) 
point or an end point of a partial orthogonal. 

The word content of each JMS entry is; 

GE635 


(36 bit word) 

J 

MU 

MD 

ISTAG 


8 

13 

13 

2 


bits 

bits 

bits 

bits 

CDC-6400/6600 

(60 bit word) 

J 

MU 

MD 

ISTAG 

24 

8 

13 

13 

2 

bits unused 

bits 

bits 

bits 

bits 


On the GE-635, the machine word length presents a practical limit to 
both the number of streamlines (J = 255) and the maximum field size 


J 

MU 

MD 

ISTAG 
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(NM - 8191) . It was felt that these limits represent a convenient maximum 
for any problem which might be encountered. Hence, the bit configuration 
for the JMS word on the CDC 6400/6600 is Identical to that on the GE635, 
resulting in an unused portion of 24 bits. The subroutines which pack and 
unpack the JMS table entries (SAVIX, GETIX, and GETBLX) are coded in the CDC 
assembly language COMPASS 1.1 in the Interest of increased computation speed. 
These routines would have to be modified if a problem required a streamline 
number J in excess of 255 or a grid point index M greater than 8191. 


6.3 STC CALCULATION STEPS, FLOW CHART, OVERLAY DESCRIPTION 

The present section outlines the sequence of calculation steps performed 
by the STC Program and describes the overlay structure on the CDC 6400/6600 
computers . 

6.3.1 STC Calculation Steps and Flow Chart 

The processing flow chart for the STC program is shown in Figure 23 and 
includes the principal subroutines, their fimction, and their associated 
output. The general calculation steps performed by the STC Program were out- 
lined in Section 5, Method of Solution. Here the calculation steps will be 
identified with the subroutine performing that step. 

The computer program has been written to analyze many types of geometries 
and flows. The user must identify the various flow bomdaries and channels 
when the program input is compiled. These sets of flow boundaries (BDY) and 
channel names (CHN) define flow regions with specified properties of tempera- 
ture, pressure, and velocity. It is the task of the program user to organize 
his problem so that each flow region is defined by the proper boimdary and 
channel names (see Figures 21 and 22) . 

The first operation in the computer program consists of reading the card 
input for the description of geometry and flow properties In each flow region 
(REDINP) , and then storing and building the various tables (BLDTAB, BLDTBS , 
BCONV) . The channel data, boundary data, and flow property data are each 
stored in their respective tables. If necessary, the boundary points are 
smoothed and local angles at each boTjndary point are calculated. 

The field point table and station table are started with the first 
unrefined grid of orthogonal and dividing streamlines. The boundary of each 
region is defined as a primary orthogonal. The dividing streamlines which 
separate the various streams are called double streamlines. 

The problem solution is started with the initial "unrefined" coarse grid 
to prevent crossing grid lines and unrealistic large streamline curvatures. 

A fully converged inner iteration is carried out to ensure reliability in the 
later stages of the calculation. 
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Figure 23. Program Flow Chart. 
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The next step in the computational procedure is to refine the very 
crude initial grid (REFINE) . A new orthogonal is placed within each region 
and a new streamline is inserted in the middle of each channel. In the 
external channel (identified by CHN EXT) , additional streamlines are placed 
close to the body. After the solution for this grid has been obtained, the 
intervals are again halved as required. 

As the calculation procedure continues, grid refinement is automatically 
provided in specified flow regions or in regions of high curvature and high 
acceleration and deceleration. The streamline and orthogonal lines which are 
added between existing lines are not required to span the field if only local 
refinement near a boundary is needed. These refinement criteria are discussed 
in Section 11 and 13. 

After refining the grid, the next step (SLC) in the solution is to 
determine the angles and curvatures of the streamlines at each grid point . 

For subsonic portions of the flow field, this is performed by fitting a piece- 
wise continuous cubic polynomial (beam) in a coordinate system which is locally 
rotated for each interval (Figure 24) . The locally rotated coordinate system 
removes the restriction that requires the slope to be small. The matching 
conditions are the angles and curvatures at each point. At the end of stream- 
lines which terminate at a flow exit boundary or extend to a flow inlet, 
the end curvature is specified. Normally the end curvature is zero, but the 
user may input a constant non-zero value of curvature. 

For grid points located in a supersonic region, the subroutine SLC 
employs backward difference formulas in keeping with the switch to the star 
with no downstream points (Figure 24). Either a 3-point parabola, a 4-point 
piecewise cubic (beam) , or a 5-point formula may be optionally selected. The 
3-point parabola is preferred. Again the coordinate system is rotated so 
that slopes in the curve fitting coordinate system are small. The end condi- 
tions for the supersonic curve fit formulas are a specified angle and zero 
curvature . 

An additional task performed by SLC is the location of stagnation points 
and the definition of the dividing streamline intersection of the boundary at 
the stagnation point. At leading edges, the dividing streamline is set per- 
pendicular to the boundary surface at the stagnation point, (Figure 25). An 
orthogonal which goes with the stagnation point is defined and the point on 
the first streamline is positioned. 

In the present version of the STC program, the angle and curvature 
calculations have been modified to improve the reliability and accuracy in the 
trailing edge regions. The special conditions enforced by these regions are 
as follows: 

a. Trailing Edge Points, Subsonic Velocity - In the previous 
version of the program, curvature at the trailing edge point had been 
taken as the boundary surface values. The curvature at this point is, 
of course, utilized in the momentum equation integration in the FLOBAL 


74 




Figure 24. Curve Fit Algorithm. 
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routine and has a strong effect on the pressure calculated at the 
trailing edge points. These pressures are utilized in iterating for the 
flow split above/below the trailing edge and obtaining "pressure closure," 
and in evaluating the trailing edge flow condition (subsonic or super- 
sonic) . A weak singularity really exists at the trailing edge (for a 
sharp wedge) , but is ignored in order to allow integration of the 
conservation equations across the entire channel and to obtain an 
approximation of the actual trailing edge pressure. 

In the new procedure, the curvature is obtained by fitting a 
parabola through the trailing edge and the two adjacent streamwise 
points. This numerical value of curvature is believed to represent 
a suitable average value for the integration of the momentum equation 
across the streamtube adjacent to the boundary. Then, in the printout 
stage of the calculation, the zero velocity, if in fact a singularity 
does exist, will be enforced. 

b. Trailing Edge Point, Supersonic Velocity - If the flow adjacent 
to the trailing edge is found to be supersonic (during the last flow 
balance Iteration), then the above calculated value is overriden by the 
boundary curvature as determined by the input boundary description. 

c. First Point Downstream of Trailing Edge, Subsonic Velocity - A 
beam is fit to the points on the dividing streamlines aft of the trailing 
edge. The end conditions on this beam are now found so that the beam 
angle bisects the trailing edge if the total pressure on the two sides of 
the trailing edge are equal. If the total pressure on one side is higher, 
than the beam angle at the trailing edge is forced to be tangent to that 
side. 


If the trailing edge has bluntness, then the wake length (stored in 
the Wake Table) is set to 2t/a where t is the "total" trailing edge 
thickness and a is the wedge angle. This allows the two dividing 
streamlines to be, in effect, circular arcs tangent to each side of the 
trailing edge and tangent to each other at the wake closing point. Also, 
beam end conditions are set to the trailing edge boundary angles. 

d. First Point Downstream of the Trailing Edge, Supersonic 
Velocity - The curvature of a supersonic point is calculated by fitting 
a beam (or parabola) to the last three points along the given streamline. 
However, in the case of the first point downstream of a trailing edge, 
a two-point curve fit is now utilized and an angle condition is enforced 
at the trailing edge point. This angle is found by calculating a 
Prandtl-Meyer expansion (or compression) to the total pressure of the 
adjacent stream if the trailing edge is sharp or to the static pressure on 
the other side of the trailing edge if the trailing edge is blunt. The 
conditions to which the flow is expanded are saved in the station table 
(as indicated above) and printed with the other output data for that stream- 
line in WRIBDY. Since the effective wake wedge angle is increased by the 
corner expansion, a new wake length calculated by the formula above is saved 
in the Wake Table. 
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PTMOVE checks the orthogonality of the grid points and moves the points 
along the spline curve as required to achieve normal intersections, CFigure 
26), between the two sets of lines. Boundary layer displacement effects are 
included at this point. Also, the normal distance, n. Is computed for each 
grid point as measured from the lower boundary of the orthogonal. 

When the near-field grid is defined, a boimdary is placed some distance 
away from the body. This boxmdary becomes the interface between the near-field 
and far-field solutions. The near-field is computed by the streamtube 
curvature method and the far-field is computed by linear small perturbation 
theory, (Figure 27). In the process of iterating, this boundary Cwhich is 
also a streamline) will float so that its shape and velocity distribution are 
matched by both the inner and outer solutions. In practice, the shape of the 
interface streamline (also referred to as the far-field boundary) is first 
assxjmed. Using the far-field equations, the velocity distribution is calcu- 
lated. These velocities computed by FARFLD, are subsequently employed in the 
near field analysis and from this comes a revised shape for the interface 
streamline. Revised velocities will then be computed in FARFLD during the 
following iteration cycle, and so forth. FARFLD is a selected boundary option 
when the far-field boundary, BDY, is called FF. A solid boundary may also 
be specified by renaming the far-field boundary and defining the coordinates. 

In the next section, subroutine ADJWF is initially called to adjust the 
flow rate at "choked" stations. For boattail analysis of nacelles, the 
internal geometry of the exhaust passage is required input to the STC program. 
Because of strea^ine curvature effects, the discharge coefficient for a 
nozzle will be somewhat less than unity. The user, however may input a flow 
rate based on unity discharge flow coefficients or, for that matter, any 
approximate value. Determination of the velocity distribution across the 
throat of a nozzle will be determined within the STC framework. The next 
three subroutines, STALOO, BRHS , and FLOBAL, develop the solution to the flow 
field equations represented by continuity. Equation 16, and radial or cross- 
stream momentum. Equation 17, (Section 5). At each station (STALOO) along a 
boundary, the flow equations are integrated along an orthogonal (FLOBAL) and 
the right hand side and the coefficient B (B = [(1 - M^) /pV^]) in Equation 20 
are determined (BRHS). In the external regions of the field, the momentxim 
equation is integrated from the far-field interface boundary to the body (or 
to the centerline or lower boundary, whichever exists). The cumulative 
streamtube flow areas are then calculated by integrating the continuity 
equation and compared with the geometric areas of the streamlines defined by 
SLC. The differences between the streamline position determined by the "flow 
balance" FLOBAL and by SLC is used as a convergence check. It also defines 
the right hand side or driving error function in the streamline correction 
equation. Equation 20. 

The next step (ADJWF) is the modification, as required, of the flow 
rates of the exhaust streams. This step is performed, if the streamline 
positions are within tolerance, by looping through STALOO, BRHS, FLOBAL, and 
ADJWF until the trailing edge departure streamline pressure compatability 
conditions ("Kutta") are satisfied. 
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Figure 27. Illustration of Far Field and Near Field Solution Domain, 


If the streamline positions are within a convergence tolerance the 
program logic loops back to REFINE. Additional streamlines and orthogonals 
are added as required by the refinement criteria. The refinement procedure 
presently built Into the program uses a criteria Involving the distance and 
velocity Increment between grid points. 

If the streamline convergence tolerance has not been met, the streamline 
correction equation is solved by defining a matrix, MCOEF, and by solving 
this matrix using one of several methods. The solution techniques include 
lAD, Implicit Alternating Direction, or LRELAX, Orthogonal or Streamline 
Block Relaxation. The user has control of the choice of solution method, 
but lAD has demonstrated the more stable and faster solutions. 

After calculating the streamline correction, 6n, the streamlines are 
adjusted (ADJSL) to their new positions and the program logic loops back to 
SLC to define new curvatures and geometric streamtube areas. This loop is 
called an inner solution. Depending on the amount of grid reflnenent, the 
inner loop may take from 2 to 20 passes to reach streamline convergence. 

As shown above, the Iterative sequence is to start with a crude grid and 
to go through subroutine SLC to ADJSL until the flow balance error is small. 
The grid Is then refined to the next level and the field is reconverged. The 
refinement /convergence process Is continued until the grid refinement criteria 
is satisfied, or alternately, until computer storage limits are reached. At 
this point, additional inner loops may be performed until the flow balance 
error is satisfactory. 

Finally, the output quantities are calculated and the results are printed 
in any of several forms. The output forms are defined in Section 8. During 
the output sequence, the boundary layer growth on the specif iced surfaces are 
determined. The boundary layer parameters <5* and d6/ds are stored in the 
boundary layer data tables for use if the STC problem is restarted. 

6.3.2 STC Overlay Description 

The STC program has been structured for execution on the CDC 6400/6600 
machines under the SCOPE 3.1 operating system. The basic OVERLAY features 
of SCOPE 3.1 have been utilized to reduce the memory requirements to those 
currently in use at the NASA Langley Research Center. Shown in Figure 28 
is the overlay structure including all subroutines and the important data 
table storage areas. As indicated in this figure, the program consists of 
a main overlay, four (4) primary overlays, and six (6) subordinate secondary 
overlays . A brief description of the processing in each overlay is given 
in the following section. In this description, the word "link" is used 
interchangeably with the word "overlay". 

Overlay (0,0) - Entry STCA 

The main overlay (0,0) contains the main program STCA, as well as general 
purpose subroutines which are called by subprograms in the subsequent primary 
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Figure 28. STC Overlay Structure. 
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and secondary overlays. The main program STCA provides control for execution 
of the STC program by loading the appropriate overlays in sequence and 
testing for completion or convergence of the problem solution. 

Also included in the main overlay is the data storage area consisting 
of the general program tables, the field tables, and the streamline table 
as discussed in Section 6.1. This storage area represents the largest block 
of working storage in the main overlay. The sizes of these various regions 
of memory are initialized by the block data subprogram USECDG. An increase 
in the number of available field points, stations or streamlines may be 
effected simply by recompilation of this block data siibprogram and the block 
data USECDM in overlay (4,0). 

Overlay (1,0) - Entry STCN 

Overlay (1,0) is the primary overlay for the program input and initial 
construction of the data tables and the unrefined calculation net. Program 
STCN serves as a driver routine to call in subordinate secondary overlays. 

These links in turn read the input and build the initial tables and calculation 
net. Subroutines used concurrently by secondary links also reside in overlay 
( 1 , 0 ). 


Overlay (1,1) - Entry STCNR 

Overlay (1,1) consists of the subroutines to read the card and/or tape 
input and perform preliminary processing. General input data are read in 
KEDINP and boundary and channel data are read in KBD and RCD respectively. 

The channel data input table is also constructed in RCD. In the case of the 
boundary input, if angles are not specified with the coordinates, the SM0TH 
routines are called to calculate angles by fitting a beam through the speci- 
fied coordinates. Also, the NACA Series 1 Cowl coordinates are stored 
internally and may be used to generate an "analytic" contour for a specified 
boundary segment. The input quantities to select these options are discussed 
in Section 8.1 and on the input sheets (Section 13.2). 

Overlay (1,2) - Entry STCN2 

The principal function of overlay (2,1) is to build the boundary table 
and the leading edge/trailing edge table using the coordinate and angle data 
from the preceding link. This procedure is acconq)lished in siobroutine BLDTAB. 
If specified by input, this overlay is recalled to initialize the matrix for 
the solution on the far field boundary (FRFDNZ) . 

Overlay (1,3) - Entry STCN3 

The bulk of initial table construction is performed in subroutine BLDTBS. 
Specifically, the orthogonal channel table, the streamline table, the station 
table, the flow adjustment table, and the wake displacement thickness table 
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are built in this routine. Subroutine BC0NV is called to build the convected 
properties table as an intermediate task. In the process of table construc- 
tion, the values of Z, R, S 2 , and V for the unrefined grid are placed in the 
field tables. 

Overlay (1,4) - Entry PRTPLT 

The principal function of PRTPLT is to provide a printer plot of the 
input boundaries and the Initial "crude" orthogonal - streamline grid. 
Orthogonals and streamlines are designated by their respective XII or XI2 
values . 


Overlay (2,0) -Entry STCB 

Overlay (2,0) is the primary link for calculation of the coefficient 
B and the RHS of the streamline correction equation using the "flow balance" 
relations. The subroutine FL0BAL and TTPT are included in this link, since 
they are used by the output routines to calctilate final velocities in the 
field prior to printing. 

Overlay (2,1) - Entry STCS 

Overlay (2,1) contains the subroutines STALOO and BRHS to loop through 
the stations and calculate the coefficients B and RHS for the matrix solution 
of the streamline adjustment equation. Included also is the subroutine 
ADJWF which adjusts nozzle flow rates before and after the flow balance 
calculation. 


Overlay (2,2) - Entry STCWl 

This overlay is loaded at the end of the problem solution to produce 
the printed output and generate a restart tape if requested. Subroutine 
WRIA prints the parameters defined by input, program storage utilization 
information, solution convergence level as well as generating the restart 
tape. Subroutines WRI0UT and WRIBDY print out the solution at the field 
points and the boundaries . If a boundary layer calculation has been specified 
for a given surface, the boundary layer output follows the normal STC 
boundary output. 


Overlay (3,0) - Entry STCXX 

Overlay (3,0) is loaded prior to the execution of the flow balance. 

The prime functions of this link are to adjust the streamline positions 
(ADJSL) , refine the grid (REFINE) , calculate streamline curvatures (SLC) , 
orthogonalize the grid (PTM0VE) , and calculate velocities on the far field 
boundary (FARFLD) . The far field routines are only called if these options 
are in effect. REFINE will not be called if grid refinement limits have been 
reached, the flow balance error is too large, or the maximum field point limit 
has been met. 
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Overlay (4,0) contains the subroutines MC9EF and lAD to set up and 
solve the maitrlx equations for the streamline adjustment. 
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SECTION 7.0 


PROGRAM NOMENCLATURE 


Communication between the subroutines in each overlay is accomplished 
by the use of labeled common. With few exceptions, the majority of labeled 
common storage areas are located in the main overlay. The size and data in 
these blocks are initialized by the block data subprograms USECDG and STCBLK. 

The principal labeled common blocks are given alphabetically in the following 
section. Labeled commons used by general purpose utility subroutines are not 
included in the tabulation. In most cases, the use of these commons is given 
in the listing of the pertinent subroutine (Voltnne II) . 

Within each block, variables are listed according to the position occupied 
in the block. In some cases, the variable namp may differ between routines and 
a typical name is given. Also, several areas are used primarily as erasable 
temporary storage and are denoted as such in the description of the labeled 
common. Inspection of the source listing should indicate the particular use 
of these areas by a given subroutine. The pertinent dimension and type infor- 
mation are included with the variable name (R = Real, I = Integer, L = Logical). 
Variables normally containing BCD data are typed as H e Hollerith, even though 
they may have real or integer names. 
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Nl I 1 Index of first point in boundary layer SW table 

NI I 1 Index of last point in boundary layer SW table 

CAPXl R 1 Equivalent flat plate distance from boimdary layer origin 

to first calculated boundary layer point. 



BLDTAl Boundary Layer Communication 

BNAMSV I 1 Current bomdary name passed to the boundary layer routines 

by WRIBDY. 
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CRXOL R 1 New orthogonal lines across a subsonic region 

CRXSS R 1 New orthogonal lines across a supersonic or mixedregion 
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Flow balance Iteration history vector 
Pressure balance iteration history vector 



JSUM I 1 Indicator used for detecting a change in the channel 

VMLBSQ R 1 Velocity squared on lower boundary 

CFKfclN Far-fleld initialization 
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SELECT Communication common between main overlay and primary 



Block Variable 

Name Names Type Dimension Description 
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SECTION 8.0 

STC PROGRAM INPUT/OUTPUT 


The following sections are concerned with the input to the STC program, 
special user instructions, and the output produced by the STC program. The 
standard system files INPUT and OUTPUT (TAPE6 = OUTPUT) are used for card 
input and printed output respectively. Additional data files designated 
as TAPEl (input) and TAPE2 (output) may be used by the programs. In general, 
these files will reside on magnetic tape. 

The standard option exercised in the STC program is to use pressure and 
temperatures in dimensionless form normalized by the free stream ambient 
pressure and temperature. When the SAB botmdary layer is chosen, however, 
pressures and temperatures must be given in dimensional form. Representative 
sets of units for input and output parameters. are given in Table I. 

8.1 PROGRAM INPUT 

Input sheets for the STC program are given in Section 13,0, along with 
special notes pertinent to the use of these sheets. Optional program input, 
not normally required, will be described in this section. Input data may be 
in the form of punched cards or a magnetic tape file (output file from a 
previous execution of the STC program). Data read from a magnetic tape file 
may be selectively over-ridden or augmented by input cards. Four (4) dis- 
tinct card input sets are read by the program and are; 


1. 

Input sheet 0 

Identification information 

2. 

Input sheet 1 

Overall input data 

3. 

Input sheet 2 

Boundary coordinates 

4. 

Input sheet 3 

Channel flow properties 


The first input set, consisting of identification information, is read 
once in a given run using fixed field format (6A10) . The remaining input sets 
consist of a header card followed by a NAMELIST input list $A. Standard 
FORTRAN IV NAMELIST (Volume II) is used to read these latter lists. Succes- 
sive cases may be run using only input sets 2 and 4. In all cases, the 
channel flow properties (input set 4) may or may not be present. 

The input parameters for each set are given in the following section. 
•Included in these descriptions are the input items appearing on the input 
sheets as well as controls for special program options and the modification 
of preset tolerances. 
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Table I. Consistent Units for STC Programs 
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8.1.1 Identification Information Input Sheet 0 

The first three (3) cards of the Input deck consist of the naine and 
address of the user and the problem Identification. 


Card No. 

Cols. 



1 

2-61 

User name (1-60 alphanumeric 

characters) 

2 

2-61 

User address or location (1-60 alphanumeric 
characters) 

3 

2-61 

Problem Identification (1-60 
characters) 

alphanumeric 


Blank cards may be substituted for Input quantities not required. 


8.1.2 Overall Input Data Input Sheet 1 




The first card of 
column 2 , the word STC 
14 and 24. 

this Input set Is a header 
starting In column 4, and a 

card consisting of a 1 In 
T or an F In both columns 

Card Column 

Description 




2 

1 - Denotes overall 

data Input 

4-6 

STC - Denotes program 

name 

14 

Input tape? (T 

or F) 



24 

Output tape? (T 

or F) 



The header card Is followed by the NAMELIST $A and 
all Input data. The NAMELIST Input Is terminated with a 

the associated over- 
$ In column 2. 

Input Parameters for General Usage 




Parameter 

Description 



Preset 

Value 

MACHO 

Mach number 




TSO 

Ambient temperature 



1.0 

PSO 

Ambient pressure 



1.0 

RG 

Gas constant 



1.0 

GAM 

Ratio of specific heats 



1.4 

RHL 

Hlghll^t radius 



1.0 

RM 

Maximum body radius 



1.0 

TTE 

Body closure tolerance 



0. 
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Parameter 


Description 


Preset 

Value 


AXI 

GR(1) 

SGR(l) 

NGR 

GZ(1) 

SGZ(l) 

NGZ 

VMGl 

VMG2 

MAXIT 

RHOC 

NODENS 


*TREF 

*MURZF 

*SCON 


Problem type T = Axlsymmetrlc T 

F = Planar (2-D) 

Table of transverse coordinates for 
grid refinement tables (up to 10 values) 

Table of transverse grid refinement 1. 

criteria (up to 10 values) 

Number of entries In the GR, SGR tables 1 

Table of axial coordinates for grid 
refinement tables (up to 10 values) 

Table of axial grid refinement 
criteria (up to 10 values) 


Number of entries In the GZ, SGZ tables 0 

Maximum Mach number Increment 0.1 

between grid points In streamwlse 
direction 

Maximum Mach number Increment 0,1 

between grid points In normal direction 

Maximum number of grid refinements 

Curvature damping factor 1.0 

Number of grid refinements for which 0 

the streamline positions are found 
by using a constant density (Based on 
total tempe,rature and total pressure) 


Reference temperature for viscosity 518.688° R 
calculation 

Reference viscosity at TREF 10 ^Ibm/ln sec 

Sutherland Constant 198.6° K 


*Used for SAB boundary layer calculation only 


Optional Input 

The following Input quantities are considered optional Input In that 
they are not normally required for execution of the STC program. These Items 
consist, In general, of Input controls for special program options and Input 
data to modify preset or Initialized constants and parameters. 
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Stagnation Properties, Optional 


Parameter 

Description 

Preset 

Valtae 

PTO 

Stagnation pressure 

1.0 

TTO 

Stagnation temperature 

1.0 

The total pressure and total temperature may be Input by specifying PTO 
and TTO normalized by the free stream static temperature and pressure. 
Alternately, the total conditions are calculated from the Mach number, MACHO 
and the ambient pressure and temperature, PSO and TSO. 


Special Controls for Supersonic Flow, Optional 


Parameter 

Description 

Preset 

Value 

SSFML 

Formula number of calculation of 
supersonic curvatures 

=1 for 3 point parabola 

= 2 for 4 point piecewise cubic (l.e. beam) 

1 


= -1 for 3 point parabolas for both 
subsonic and supersonic points 


SSEF 

Supersonic entering flow along the 
upstream boundaries? T or F 

F 

SSEANG 

Entering flow angle for SSEF = T (degrees) 

0. 

SSDF 

Supersonic flow downstream of a choked 
station? T or F 

F 

SSFEND 

Stxpersonlc beam downstream end 
condition (SSFML =2) 

.75 


0 . = Parabola 

1. = Cubic 


SSFMDl 

Supersonic beam upstream end 
condition (SSFML =2) 

.75 


0. = Parabola 

1. = Cubic 


TSIC 

Number of points to be tjsed for 
transonic interpolation of curvature 

2. 

KHOW 

Flow difference damping factor 

1. 

RHOWSS 

Supersonic - flow difference damping factor 

1. 

BHOCSS 

Supersonic - curvature daiq>lng factor 

1. 
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Boimdair Conditions for Streamline End Points, Optional 

Preset 


Parameter 

Description 

Value 

NBCIN(l) 

NBCIN(l) E Upstream, NBCIN(2) = Downstream 
NBCIN = 1 Angle Specified 

NBCIN - 2 Curvature specified 

2 

ACF(l) 

ACF(l) E Upstream, ACF(2) E Downstream; 
Angle or curvature 

0 . 

Matrix 

Solution Controls, Optional 


Parameter 

Description 

Preset 

Value 


I ADM Sweep parameter 0 

0 = lAD 

1 = Orthogonal block relaxation 
-1 = Streamline block relaxation 

RHOBAS Base acceleration factor 0.5 

parameter, Pg 

RHOAMP Half amplitude of the acceleration 0.5 

factor, 

Note: 

^ „ , 2 fniT ' 

p = P-, + 2p. sin — 

MAXSWP Maximum number of sweeps 200 

TOLRL Relative tolerance for matrix solution; .001 

1 ADS 2 I < T0LRL*DS2 

' ' — Tnav _ 


Flow Balance Solution and Flow Adjustment Controls, Optional 


Parameter 

Description 

Preset 

Value 

T0LES2 

Relative tolerance on maximum point 
movement predicted by the flow balance 

0.001 

TOLINR 

Inner iteration tolerance on streamline 
movement 

0.05 

NINNER(l) 

Number of inner iterations (without 
grid refinement) . Specify up to MAXIT 
values 

16*10 
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Preset 

Parameter Description - Value 

CNVF(l) Fractional percentage of total point 10*1.0 

movement to be used at a given refinement 
level. Specify up to MAXIT values 

TOLWF Tolerance on the fractional flow 0.001 

adjustment needed to meet the trailing 
edge pressure closure condition 

Grid Refinement Parameters, Optional 

The following Input Items may be ttsed to control the length of ortho- 
gonal lines and streamlines as the grid Is refined, the use of these 
parameters Is discussed In the following section on "special user Instructions 


Preset 

Parameter Description Value 


CRXSL Used for extending new streamlines .375 

CRXOL Used for extending new orthogonal .375 

lines across a subsonic region 

CRXSS Used for extending new orthogonal .125 

lines across a stipersonlc or mixed 
region 

CRXE Used for extending new orthogonal 0 

lines which cross an expanding 
flow sonic line 

CRXC Used for extending new orthogonal 0 

lines which cross a supersonic to 
subsonic compression line 


Optional Print and Control Options 

Optional diagnostic print and special logic controls are located In 
common /CPRINT/. Normally these variables are set to 0.. As Indicated In 
the following section, the print controls will provide output from a number 
of routines. The setting, pertinent subroutine, and resulting printout are 
Indicated. Several variables control logic flow only. 

Parameter Description 

PRTES2, PDD(l) Subroutine BRHS 

=2 Print B,RHS,DS2,Z,R,PHI,CURV, and ES2 
for the two stations with maximum ES2. (DS2 
Is the value con?)uted during the previous 
Iteration.) 

>2 Print the above information for 

PDUM(8) PRTES2. (B,RHS and DS2 are 

the values computed during the previous 
iteration.) 


107 



PDD(3) 


Subroutine lAD 


PDD(6) 

PDUM(l) 


PDUM(2) 


PDUM(3) 


#0 Print matrix solution for DS2 at 
end of each sweep. This Indicator 
Is turned on Internally If 
NSWP > MAXSWP - - 2 

Sxibroutlne BRHS 

=2 Suppress area scaling for ES2 
calculation 

Siibroutlne SLC 

>1 Print 2ffi,RB,ANGD,CURVD during 

Iteration for stagnation point location 

=1 Print results of curvature calculation 

=2 Print results of supersonic curvature 
calculation 

=4 Print curvature results If point Is a 
stagnation point or the termination of 
a partial orthogonal 

Subroutine SLC (Sharp T.E, Logic) 

=0 Treat T.E. singularity by numerical 
approximation. Average angle and 
curvature are obtained by a 3-point 
curve-fit and these valiies are printed 
in the tabulated WRIBDY output. 
Supercritical flow on the other side 
of the T.E. will be expanded to the 
confuted (static) pressure, like the 
blunt T.E. case 

(=1) Use the same logic except that a 

Prandtl Meyer fan, if one exists, is 
expanded only to the total pressure. 

On the concave comer side, the actual 
flow angles and stagnation pressure 
are printed in WRIBDY Instead of the 
nximerically computed values . 

Subroutine SLC (Wake decay rate) 

=0 Set wake angle to physical wedge angle 
minus Prandtl Meyer expansion angle. 

0 < PDUM(3) <1 is a danq>lng factor; 
viz, 0 = no damping. 

=1 Use physical wedge angle. 
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-PDUM(4) SubrOutlae SLC (curvature formula dovmstream of T.E.) 

”0 Use the "End Internal Parabola" (FEND •= 0) for end 
condition on the trailing streamline. Also, use 
3-polnt backward curvature formula for 1-st point 
downstream of T.E., If M > 1. 

(=1) Use T.E. angle plus Prandtl Meyer expansion as 

an end condition for a 2-polnt curvature for^la If M>1. 

=2 Use downstream T.E. singularity angle for the beam curvature 
formula (M<1) . 

(Set PDUM(2)=1 for compatible printout) 

PDUM(5) Subroutine MCOEF 

=1 Print Influence coefficients, G, relating 
streamline point movement to the negative 
of curvature 

=3 Print above Information for partial 
orthogonal points only 

=4,5 Print above Information for 1st and 2nd 
points upstream of leading edge 

=4,5 Set FEND(2)=0 for leading edge stagnation point 
(FEND (2) =1 standard) 

PpUM(6) Subroutine ADJWF 

>0 Print results of flow adjustment Iteration 

PDUM(8-9) Subroutine FLOBAL 

Print flow balance related data If 

PDUM(8) PDUM(9) and 

PDUM(8)*PDUM(9) # 0 

PDUM(IO) Subroutine BUILDT 

=1 for printout of tables described In 

Section 4.1 after they are first constructed 

=2 for printout of all tables at end of 
solution 

PDUM(ll) Subroutine FARFLD 

#0 Print coordinates, slope, velocity on 
the far field boundary 
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Parameter 


Description 


PDUM(12) 

PDUM(13) 

PDUM(17) 

PDUM(18) 

PDUM(20) 

PD^(19) 

PRPRN 


Subroutine MCOFF 

00 Use spline Influence coefficients even 
If SSFML = -1 

Subroutine FARPLD 

=0 Use "linear" Interpolation on the far field 
boundary 

00 Use "least squares parabolic" Interpolation 
on the far field boiandary 

Subroutine BDYPTM 

=0 Print separated boundary layer message 

9^0 Delete print of separated boundary layer message 

Subroutine FLOBAL 

=1 for calculation of normal shock total 
pressure loss for transonic compressions 

Subroutine FRFDNZ 

00 Print matrices used In far field solution 

Subroutine SLC 

=2 Omit stagnation point Iteration on first 
time through SLC 

Subroutine WRIOUT 

0-1 Print field information at each station 


Arbitrary Pressure and "Free" Boundary Input, Optional 
Parameter Description 


Preset 

Value 


NZP 

Nuinber of entries in tables of ZP vs PPS 
for pressure and free boundaries. NZP 

NZPl 

Number of entries which apply to first 
boundary if two are present 

ZP(1) 

Table of axial distance 

PPS(l) 

Table of pressures 

PSPIV 

=0 PPS entries are pressures 
=1 PPS entries are velocities 
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The special boundary names which trigger the pressure boundary use are 
described In a following section. When two boundaries are used, entries for 
the second boundary follow the entries for the first boundary. When a 
section of the flow field Is bounded by two pressure or free boundaries, the 
end boundary angle must be specified. 


ACF(2) = _, NBCIN(2) = 1, 


Blockage/Lamina Thickness Input, Optional 


Parameter 


Description 


NTHKX 

THKX 


Number of entries In the THKX table 
2 < NTHKX £ 25 (axial) 

X or Z axial coordinate 


NTHKY 

THKY 

THIK2D(J,I) 


Number of entries In the THKY table 
1 £ NTHKY < 25 (vertical) 

Y or R vertical coordinate 

Fraction of unblocked circumference or 
lamina thickness THIK2D <250 


Preset 

Value 


In the THIK2D table, the vertical variations must be listed first; viz, 

ITHIK2D (J,I), J - 1, NTHKYJ, 1^1, NTHKX 

The THIK2D table Is not extrapolated. Outside the range of THKX or THKY, 
the end values in the THIK2D table will be used. 


Boundary Layer/STC Restart Cases 


Parameter 


Description 


Preset 

Value 


INPBLR Numbers of Input boundaries for which 0 

a boundary layer calculation is to be 
performed 


When INPBLR Is Input, INPBLR fixed format cards must follow the 

$A $ overall Input data namelist set. The Information on these cards 

consists of the following Items: 


Card Column 


Description 


2-7 

12-21 


Boundary name, 1-6 alphanumeric characters 
Equivalent flat plate distance from boundary 
layer origin to the first calculated boundary 
layer point (F10.6 format) 
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The second input item may also be specified as CAPXl input list (see 
next section). The boundary data list may be read only once ; hence, it is 
normally necessary to specify a boundary layer calculation at the beginning 
of the problem. In many cases, it is desirable to obtain a fully converged 
STC solution before introducing botmdary layers into the calcxilation. The 
above input allows the introduction of bovmdary layers at any stage of the 
solution by using a restart procedure and following the restart $A Namelist 
with the fixed format boundary layer input. 

8.1.3 Boundary Coordinates Input Sheet 2 

The first card of this input set is a header card containing a 2 in 
column 2, the word BDY, starting in colisnh 4, the boundary name in colxjmn 14, 
and the channel name in column 24. 


Three special boimdary options are provided: far-field, pressure, and 
free. The far-field option may be used for an external flow upper boimdary 
and may be Invoked by naming the boundary FF. This name causes the 
numerical flow solution to be matched to a small perturbation analytical 
solution in the region from the "far-field" bomdary to infinity. 

An arbitrary static pressure may be specified along one or two boundaries 
by using the boundary names PEESl and/or PRES2. The boundary pressure 
specification is described in Section 8.1.2. 

A free boundary is a constant pressure boundary which is downstream of 
a fixed boundary and has the same pressure level as the pressure at the last 
point on the fixed boundary. The special names for this type of boundary 
are FREEl and FREE2. Further Infor^tion on this type of botmdary is given 
in the input sheet notes CSs-ction 13l. 


Card Column 


Description 


2 

4 

14 

24 


2 - Denotes boundary input 

BDY - Denotes input type 

Boundary name, 1-6 alphanumeric characters 

Channel name, 1-6 Alphanumeric characters 


The header card is followed by the NAMELIST $A and the associated input 
for the specified boundary. The NAMELIST is terminated with a $ in column 2. 


Parameter 

Description 

UPPER 

Boundary position 



UPPER = T 

Upper "bomdary 


UPPER = F 

Lower Boundary 

ZRONLY 

Geometry indicator 



ZRONLY = T 

No surface angle input 
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ZRONLY = F 

Surface angle input 


Two options are available for the Input of the boundary geometry. Either 
the coordinates and the surface angle (measured from the positive x-axis) may 
be Input, or the coordinates alone may be Input. In either case the points 
must be input accurately. The first option (ZRONLY = F) is preferred. 

With the second option (ZRONLY = T) , a beam is fit to the input points to 
determine the angles and, in the process, the beam fit angles and curvatures 
are printed. To determine whether a suitably smooth curve has been fitted 
to the points, the user should examine the curvatures and make sure they 
are reasonable. In general, the points should not be closely spaced except 
as required in regions of hi^ curvature. In these regions, the angle 
change between points should be less than 25 degrees . 

The NACA Series 1 cowl coordinates are stored internally, and may be 
selected by specifying ZRONLY = T and B(l) = 991, 1, followed by the high- 
light coordinates and the coordinates of the maximum diameter. , If the cowl 
is to be extended beyond the end of the Series 1 contour, the "maximum 
diameter" coordinate is repeated and then other coordinates are listed. 

Boimdary coordinates are normally input in tabtilar form using the B 
block as specified on input sheet 2. All points must be listed in the stream- 
wise direction for each boundary. Points at sharp comers must be listed 
twice, once for each angle which exists at the point. Normally, pressure 
and Mach number distributions will be printed at each orthogonal intersection 
with the boundary. Orthogonals may be forced to coincide with boundary points 
by listing the point twice in the input and setting DBLPTS =0. (see optional 
boundary input) . 

The column names in the B array are Z (or X) , R (or Y) , and ANGD, 

Parameter Description 

B(l) Input block for boimdary data 

column 1 Z or X coordinate 
column 2 R or Y coordinate 
column 3 Angle or slope of surface 

If desired, the data in the B array may be input in "free form" using the 
symbolic nan^s associated with the appropriate columns of the B block. 


Parameter 


Description 


R 

Y 

Z 

X 

ANGD 


Normal (radial) coordinate (axisymmetrlc) 
Normal (vertical) coordinate (Planar) 
Axial coordinate (asix 3 rmmetric) 

Axial coordinate (planar) 

Angle or slope of surface 


Boundary data input via the B block will override corresponding data 
input in the "free form". 
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Boundary Layer Specification 


The preset program option is to not calculate a boiindary layer on a given 
surface. 

The necessary input to specify a boundary layer calculation is as follows: 

Preset 


Parameter 

Description 

Value 

BL 

F - No bomdary layer 
T - Bomdary layer 

F 

CAPXl 

Equivalent flat plate distance from bomdary 
layer origin to the first calctilated 
bomdary layer point. 



Optional Input 

Optional input parameters, not shown on input sheet 2, are available to 
apply linear transformations to the input coordinates. Also, input control 
may be specified to force orthogonals to be placed at selected boundary 
input points . This is accomplished by repeating the coordinates in the 
bomdary input and specifying DBLPTS = 0. Normally, DBLPTS = .01 forces 
removal of extra orthogonal stations where the angle discontinuity is less 
than .01 degrees. 

Preset 

Parameter Description Value 


DBLPTS 

Double point tolerance. Double points 
will be deleted if the angle discontinuity 
is < DBLPTS (degrees) 

ROTATE 

Angular rotation in degrees 

ZPIVOT 

RPIVOT 

Coordinates of the pivot point of 
rotation 

SCALE 

Multiplicative constant to be applied to 
the coordinate data 

ZTRANS 

Translation incren«nt in axial direction 

RTRANS 

Translation increment in radial/vertical 
direction 


8.1.4 Channel Flow Properties Input Sheet 3 

The first card of this input set is a header card containing a 3 in 
column 2, the word CRN starting in column 4, and the channel name starting 
in column 14. 
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Card Column 


Description 


3 3 - Denotes channel Input 

4 CHN - Denotes Input type . 

14 Channel name, 1-6 alphanumeric characters 


Channel flow properties are specified in the following NAMELIST 
is terminated with a $ . 

A, which 

Parameter 

Description 

Preset 

Value 

GAM 

Ratio of specific heats 

1.4 

RG 

Gas constant 

1.0 

TTO 

Total temperature normalized by 
free stream static temperature 


PTO 

Total pressure normalized by free 
stream static pressure 


MACHO 

Mach number 


TSO 

Static tenq)erature normalized by 
free stream static temperature 

1.0 

PSO 

Static pressure normalized by free 
stream static pressure 

1.0 


Specify either TTO and PTO (If known) or MACHO, ISO, and PSO If the Mach 
number Is known. When these options are not iised, the second option Is 
employed with the free stream values of PTO and TTO computed from the MACHO 
supplied li? the overall Input data (sheet-1) . 

Preset 

Parameter Description Value 

AO Nondlmenslonal flow area normalized by 

the product of ir and the square of the 
highlight dimension. For plane cases 
Input AY /AY„- . For Inlets AO = mass flow 
ratio. ° 

The flow Is computed using the total properties as determined from the 
supplied Mach number, MACHO, and the flow area AO. If input data are not 
supplied for a given channel, the reference properties on Input sheet-1 
will be employed with the area calculated at the entrance station. 

Individual sets of bomdary and channel data need not be input in any 
particular order. For exarq>le, boundary Input for the upper and lower sur- 
faces bounding a given channel (input sheet 2) may be followed by the input 
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for the pertinent channel (Input sheet 3) . While arbitrary placement of 
these input sets may be used, it is well to develop some standardized con- 
ventions for the sequential input of boundaries and channels. This will 
be discussed further in the following section on special user instructions. 


8.2 PROGRAM OUTPUT 


The output from the STC program may be logically divided into the 
following eight sections . 

1. Card input and preliminary printout 

2 . Input and calculated boundary coordinates and angles 

3. Printer plot of initial XII, XI2 grid. 

A. Solution history 

5 . General input and output data 

6 . Flow field data along orthogonal lines 

7. Calculated flow data along field boundaries and final channel 
momentum balances. 

8. Boundary layer data 

The above sections appear sequentially in the output except when 
PRPRN = -1. In this case, the print of the flow field data along orthogonal 
lines is eliminated and the calculated boimdary information and boundary layer 
output follow the general input/output section. 

8.2.1 Card Input and Preliminary Printout 

The initial section of output consists of a card image print of the 
problem input and a designation of the tape Input/output file selections; 
viz, TAPIN = T/^ and TAPOT = T/F. Upon completion of the card image print, 
the file TAPES is rewound to its original position. 

8.2.2 Input & Calculated Boundary Coordinates & Angles 

Boundary coordinates are processed and printed as they are read. The 
columnar output displayed for each boundary consists of the input coordinates 
(X,Z), (Y,R) and the input or calculated body slope (ANGD) in degrees. 
Preceding these items is a bold heading specifying BOUNDARY COORDINATES, the 
boundary name (BDY = name) , the adjacent channel (CHN = name) , whether the 
surface is an upper (UPPER = T) or a lower (UPPER = F) boxmdary, and boimdary 
layer specification (BL = T,F) . When the ZRONLY = T option is in effect, 
intermediate printout will be produced defining the results of the beam curve 
fit to the input coordinates. Since no smoothing is applied to these 
coordinates, the only meaningful output is the (X,Z) (Y,R) and the calculated 
angles (ANGD) . Break points in boxmdarles (double points) enable the curve 
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fit of a boundary in segments. The consolidated output, consisting of the 
collated boundary data for all segments, follows the above intermediate 
printout. The NACA Series 1 cowl coordinates (ZRONLY = T) appear as 40 points 
in a form similar to that described previotosly. 

Normally, a plot of the XII, XI2 initial grid (described in Section 
9.2.3) is printed as the next section of output. When far-field boxmdary 
conditions are applied, however, a comment is inserted on the output designat- 
ing the orlgineil transverse location and streamwise extent of the far field 
boundary. Prior to calculation, the upstream and downstream limits of the 
far field streamline are extended to insure that the velocities are well 
behaved at the streamwise ends of the STC integration region. The coordinates 
of the resulting far field boundary are printed below the above comment. 

8.2.3 Printer Plot of Initial XII, XI2 Grid 

The STC program uses an auxiliary nvnnbering system (Ci»C2) for orthogonal 
lines and streamlines. The coordinate = XI and ^2 = ^12 apply to the 
orthogonal lines and streamlines respectively. A printer plot of the initial 
XII, XI2 grid, along with the nominal bomdary shape appears as the next 
section of output. 


8.2.4 Solution History 

The solution history output provides a convenient summary of the problem 
history as the solution proceeds through major grid refinements and inner 
iterative improvements of the flow balance. Included in this section are the 
parameters defining the grid refinement, the inner iterations, the matrix 
solution, the flow balance error, and the flow adjustment parameters for the 
"Kutta" iteration at trailing edge points. The pertinent variables in their 
literal order of appearance in the solution history list are: 

Variable Description 

NREFIN Grid refinement level 

Number of field points 

Inner (flow balance) iteration level 

Number of imbedded supersonic points 

Number of iterative sweeps for solution of 
the matrix equations. 

= 0 when INRCTR = 0. 

MAX-DS2 Maximum streamline correction predicted 

by the matrix solution for DS2. 

= *0.000000 when INRCTR = 0. 

MAX-ES2 Maximum ES2 predicted by flow balance 

LIM-ES2 Limiting streamline ES2 for satisfactory 

solution of the flow balance. 


GRID PTS 
INRCTR 
NSSPTS 
NSWEEPS 
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Variable 


Description 


*Tralllng 

Edge-XI2 

*Flow 

Rate 

*Fractional 
Flow Error 


* This output is blank if there are no trailing edges 
8.2.5 General Input and Output Data 

The output in this section consists of general input and a. summary of 
selected output parameters. All of these items have been defined previously 
in Sections 8,1 and 7.0, hence only the general subdivisions of this printout 
will be listed. 


XI2 coordinate at trailing edge 

Variable channel flow rate 

Fraction flow error during flow adjustment 
Iteration for satisfaction of the "Kutta" 
condition 


a. General Input 

b. Streamline End Conditions 

c. Supersonic Curvature 
Parameters 

d. Subsonic/Supersonic 
Branch Selection 

e. Grid Size Criteria 

f. Memory Utilization 


Section 8.1 
Section 7.0 
Section 7.0 

Section 7.0 

Section 8.1 


The memory usage in terms of field points, STC table storage, and 
number of streamlines is con^iared with the maximum available values 
of these parameters. 


g. Convergence Data Section 8,1 

h. Special Boundary Options 

i. Matrix Solution Section 7.0 

j. Contents of Channel Input Table Section 7.0 

k. Channel Flow Rates, Pressures 

and Temperatures 


Included in this output are the channel flow rates and the pressures and 
temperatures in each individual channel. Note that the flow rate in a given 
channel may be adjusted to satisfy the Kutta condition at a trailing edge. 
When dimensional properties are input, the columns denoted PT/PSO and TT/TSO 
have units compatible with the dimensions of the input gas constant RG. 
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8.2.6 Flow Field Data Along Ortho Ronal Lines 

Flow field data on orthogonal lines are printed when PKPKN // -1. The 
first portion of output on each page consists of the program name (STSEAMTUBE 
CURVATURE PROGRAM) and the problem identification. Following the problem 
identification is a heading with the constant station value = XII along with 
the channels through which the given orthogonal passes . Primary (Initial) 
orthogonals are flagged with ** after At the extreme right of this line 

of output is the flow type; viz., SUB = subsonic, SUP - supersonic, or CHOKE = 
choked. Pertinent data from the field arrays are printed in column format. 
Stagnation points do not appear in the flow field output. 


Variable 


Description 


X12 

STEM FNCT 

X, Z 

Y, R 
PHI 
CURV 
PS/PO 

PS/PT 

TS/TT 

CP 


Streamline numbering coordinate 5 2 

Dimensionless stream ftmctlon; fraction of flow 
in channel 

Streanwise coordinate 

Transverse or cross stream coordinate 

Flow angle, deg. 

Streamline curvature 

Ratio of local static pressure to reference 
static pressure 

Local static to total pressure ratio 
Local static to total terq>erature ratio 
Pressure coefficient 


Compressible 


Constant density 

(y = 0) 

0 MACHA < . 1 

MACH 
AREA 
PTQPTO 
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Mach number 
Flow area 

Channel total pressure ratio 


P - P - 
s “ 

|yP M^ 

Z 00 oo 


p - P 
s " 


PM 

OO 00 



For primary orthogonals, the axial and normal momentum £l\ixes are 
printed following the field data. These items are: 




Axial Momentum Flux = fV cos (j» dw + l(P-P ) cos 4> dA 

s 


Normal Momentum Flux 




sin (|) dw + 


/< 

/'■ 


P ) sin (J) dA 
s 


8.2.7 Calculated Flow Data Along Field Bouhdaries-Final Channel Momentian 
Balances 


The final section of standard STC output includes the flow parameters on 
the upper and lower streamlines bounding each channel. Normally the bomdary 
data follows the field data. In the case of PRPRN = -1, however, the boundary 
information is printed immediately after the general Input/output section. 

Three lines of identification information head the bomdary data, and 
consist of the program title (STREAMTUBE CURVATURE PROGRAM) , the specific case 
Identification, and a designation as to whether the boundary is an upper or 
lower portion of the given channel. The streamline coordinate (XI2) is ’ 
given following the channel name. Boundary flow parameters are printed at 
orthogonal intersection and consist of the following items: 

Description 

Orthogonal (station) numbering coordinate (C^) 

Distance along boundary 
Axial coordinate 
Vertical or radial coordinate 
Flow angle or surface angle 
Streamline or surface curvature 
Local static to reference pressure ratio 
Pressure coefficient 

Cp - 


Variable 

XII 

SIW 

XW, ZW 

YW, RW 

ANGW 

CURVW 

PS/PO 

CP 


C 


P 


= R8(P-P,)/P,m/ 


(Y = 0, constant 
density) 



(“o 


< . 1 ) 
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Variable 


Description 

PS/PT Local static to total pressure ratio 

MACH Local Mach number 

PT/PO Ratio of local total to reference static 

pressure 

The cumulative drag/ thrust Is printed In the column labeled CDPI. The 
force Is normalized by the free stream dynamic pressure (q) and the maximum 
area based on EM and Is given by: 


CDPI = 


("'-''ref’ 


dA 


<1 A 


max 


The adjacent column (A-AMAX) /AMAX Is the projected area normalized by the 
maximum area. Following the column printout Is the ratio of boundary total 
temperature to ambient total ten^erature . When the boundary Is the approach 
stagnation streamline, the Inlet additive drag Is printed below the total 
temperature ratio. This force Is also normalized by q and AMAX and Is printed 
for both the upper and lower stagnation streamline. 

Momentum Balance 


The STC program evaluates the thrust/drag on each boundary surface 
bounding a given channel and then verifies these forces by performing an 
"overall" momentum balance for each of the fluid streams. The Integral 
momentum balance output follows the boundary output for a given channel. This 
Information consists of the entering axial momentum, the Integrated pressure 
forces on the boundaries, and the leaving axial momentum. The discrepancy 
between the leaving momentum and the sinn of the Inlet momentimi and the 
pressure- area forces represents the net error In the calculation. This error 
may be attributed to Inaccuracies In the coiiq>uted pressure distributions or, 
perhaps, to Insufficient refinement of the calculation grid for adequate 
resolution. 

{ 

8.2.8 Boundary Layer Data 

The standard output from the SAB portion of STC consists of the boundary 
layer parameters at each orthogonal Intersection of the boundary. The 
Initial output consists of a bold heading specifying BOUNDARY LAYER. This is 
followed by the botmdary layer parameters. 
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Variable 


Description 


XW 

THETA 


DSTAR 


Axial coordinate 
Momentum thickness 



Displacement thickness 


DELTA 

REX 

CAPX 


CF 


SW 

DSTR 




dm 


Boundary layer thickness 6 

- p V s 
0 0 

Local Reynolds Number R^^ = — 


Equivalent flat plate distance along surface 
S 


Pr^Jl 


Pr^ds 


where P = 


M 


1 + ,2M^ 


a = 0 Planar flow 
a = 1 Axisymmetric flow 

Skin friction coefficient 


C. = 


w 


1/2 P V 

e e 


T = Shear stress at solid surface 


w 


Distance along surface 
"Smoothed" displacement thickness 6 
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Variable 


Description . 

DDSTR Derivative of DSTR (dfi /ds) used for correcting the 

local flow angle 

SEP Separation flag. Normally -blank, appears as SEP if the 

boundary layer has separated 

F Stratford separation parameter 

The Stratford separation parameter is defined as follows: 
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SECTION 9.0 


PROGRAM MESSAGES AND ERROR CODES 


9.1 PARAMETERS DESCRIBING SOLUTION STATUS 


The solution history output, described briefly in Section 8.2.3, contains 
a summary of how the solution is proceeding. By reviewing this solution 
status, the level of grid refinement and the degree of solution convergence 
can be readily determined. In brief, the solution history is like a "hospital 
patient's history" in that the general health of the calculation procedure 
is charted here. (See Section 12 - Sample Cases). 

The grid refinement level is specified by the number of times the grid 
has been refined (NREFIN) and by the number of GRID PTS in the flow field. 

The number of grid refinements is controlled by MAXIT. The number of points 
in the flow field is dependent on the grid refinement criteria specified and 
MAXIT. 

The number of inner iterations at any grid refinement level is shown by 
INRCTR. These are controlled by built-in tolerances on the flow balance, 
but the internal program control can be overridden by inputting NINNER(m) 
where m refers to the level of grid refinement. 

The system of solution tolerances is shown schematically in Figure 29. 

The iteration tolerances are controlled by TOLINR, T0LES2, and TOLWF which 
may be input. The system starts with the maximum streamline movement 
(MAX-ES2) demanded by the "flow balance" in FLOBAL and compares to a specified 
tolerance, the "flow balance" is satisfied. At this point, the flow rates 
of variable channels are adjusted to satisfy the "Kutta" condition at all 
trailing edges by iteratively calling for a "flow balance" and a "flow adjust- 
ment" (FLOBAL and ADJWF) . The summary of the flow adjustment procedure is 
shown under the heading "KUTTA ITERATION" in the solution history output. When 
no trailing edges are present, the flow adjustment procedures are bypassed, 
and the solution immediately calls for additional grid refinement (if the limit 
has not been reached. If the MAX-ES2 is not within limits, the inner Iterations 
continue until the specified tolerance, LIM-ES2, is met ''(or the specified 
NINNER limit is reached). The specified tolerance, LIM-ES2, is controlled 
by TOLINR and the current grid size. Thus the "flow balance" tolerance 
becomes smaller as the grid refinement proceeds. 
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Note: Note; 

LIM-ES2 = TOLINR* (Current Grid Size) HAXIT, NINNER, TOLINR, TOLES2, and 

CLEN = Characteristic Grid Interval Length; TOLWF are Input Items 

Average of SGR's & SGZ's 


Figure 29. Iteration Tolerance. 

If the grid refinement (MAXIT) criterion is satisfied or the field point 
limit is reached, a second tolerance (T0LES2) defines a final series of inner 
iterations. This permits a looser tolerance (TOLINR) during the flow field 
grid development and then a final tight tolerance (T0LES2) for actual flow 
field definition. The summary of this iteration tolerance is shown as NREFIN, 
GRID PTS, LIM-ES2, MAX-ES2, Z, and R. The latter two parameters are the field 
coordinates of the maximum streamline movement point. A comparison of 
MAX-ES2 versus LIM-ES2 (See Section 12 - Sample Cases) normally shows decreas- 
ing values of MAX-ES2 as the solution proceeds. Any erratic behavior in 
MAX-ES2 indicates convergence problems. 

During an inner iteration, the streamline displacement (or correction) 
equation is solved to adjust the streamlines in the flow field. The results 
of the matrix solution are shown as a maximum streamline adjustment (MAX-DS2) 
which is actually made and the number of sweeps (NSWEEPS) throu^ the matrix 
solution (lAD or LRELAX) to reach the final solution. 
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The final parameter in the solution history output is NSSPTS, the number 
of imbedded supersonic points in the flow field. This item, in addition to 
the erratic behavior of MAX-ES2 may be indicative of solution convergence 
problems . 

Under certain output options, the diagnostic printout Includes a line of 
solution history at the end. Since the history is not summarized neatly on 
one page, care must be exercised in seeking out the solution history inter- 
spersed with the diagnostic output. 

9.2 STC ERROR PROCESSING 


In general, the occurrence of an error during processing will cause 
termination of the problem. Each of the main links contain an error print 
routine which will print the STC tables and the field tables to assist in 
diagnosing the problem. When the error occurs in overlays subordinate to 
(0,0), an illegal computed GO TO statement is executed. This will cause a 
MODE 4 abort and a printout of the sequence of subroutine calls preceding the 
error. 


9,2.1 Description of Table Print Format 

In general, the information stored in the STC tables consists of mixed 
BCD, floating point, integer, and logical data. A special routine (TABPRT) 
is used to print the arrays which may have a variable data format. The 
following conventions have been established for the output format on the CDC 
machines : 


Data Type 


Print Format 


Integer 

Floating Point 
BCD 

Logical 


Integer (112) Format 

E or F Format, depending on the 
magnitude of the number. 

Hollerith characters A6 Format 

TRUE. Variables are printed as an 
integer -0 

FALSE, Variables are printed as ah 
integer 0 


Note that a floating point 0, an Integer 0, and a logical. FALSE, are all 
printed as an Integer 0 since the data type, in this Instance, is indisting- 
uishable. The program listings (Volume II) may be consulted to determine 
the true type of the data should this situation arise, A junk word BITS = 
10^5 Is utilized in the STC program for testing purposes. This quantity is 
printed using OCTAL format (04) and appears as 0000. 
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Common/TRDUBL/contalns two logical variables ERR and INERR which may be 
set to. TRUE, during processing. In this sltutatlon, the job Is terminated 
normally In the main program. The diagnostic output Is preceded by the 
comment ERR = T» ERRLOC = ni, LENTRY = n 2 * In this comment, n^ and n 2 are 
Integers denoting the location of the last overlay call In the MAIN program 
of OVERLAY (0,0) and the current LENTRY setting (see listings) . When an error 
Is detected In overlays subordinate to (0,0), a mode 4 error Is Initiated to 
terminate the program with a subroutine back trace following selected diagnostic 
output. 

Since the data tables used In each overlay differ, the error print 
may differ slightly between links. Normally, common printed output consists 
of the general labeled common /ALLCOM/, the STC tables (BDYTAB, STATAB, etc.), 
and the field tables which Include R, Z, VM, etc. The field tables are 
printed in tabular format by point number (M) for ease in locating computation 
prob lems . 

In the STC execution overlays (2,n) and (3,n) , an additional block of 
output Is printed, consisting of the information currently in the /ERASE2/ 
common block, used as temporary storage by the major STC subroutines. These 
data are headed by the subroutlne(s) currently using the temporary storage 
and the variable names currently applicable for common /ERASE2/. 


Subroutines! BRHS. FLOBAL. WRIBDY. WRIOUT 


Data are printed in tabular format along the current orthogonal line 
(subscript K = 1,NK) 


Parameter 



AREA 

AREAO 

DISP 

PT 

LAMBDA 

RHO 

SQRTW 

TS 

TT 

VMSQ 


Flow area calculated during FLOBAL Iteration 

Initial FLOBAL streamtube flow areas calculated from S2 
tables. AREAO(NK) Is the passage flow area. 

DISP(K) is the displacement thickness of the wake between 
streamlines K,K+1. DISP(K) = -1 for a zero thickness 
wake . 

Total pressure 

Fraction of unblocked circumference (axisymmetric) or 
lamina thickness (planar) . 

Density 

Square root of WKQKP 
Static temperature 
Total temperature 
Velocity squared 
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WKQKP 


Ratio of velocity squared on streamline K to value on 
streamline K+1 due to curvature and total tenq)erature 
change 

WQA Flow per unit area, pV 

WSTA Cumulative flow rate as obtained from streamline table, 

SLTAB. 

RG Value of the perfect gas constant for each streamline 

C2CP Two times the specific heat at constant pressure 

FGR Value of 1/(y”1) for each streamline. 

Subroutine; SLC 

Data are printed In tabular format along the current streamline 
(subscript 1=1, NI) 

Parameter Description 

CURSS Array for calculation of streamline curvature at super- 

sonic points 

QV Stagnation point Iteration control vector 

RB Radial distance, ordinate, for points along streamline 

ZB Axial distance, abclssa, for points along streamline 

ANG Flow angle for points along streamline 

CURVB Curvature for points along streamline 

SIB Distance along streamline 

BI Subsonic - supersonic Indicator for points along stream- 

line 

J2D0NE Indicator for streamlines which have been curve fit. 

Used for Indicating delayed Interpolated data at partial 
orthogonal termination points, 

MSV Field point Index along streamline 

Sub routine ; PTMOVE 

values along the control streamline (I = l.NIC) 

Parameter ' Description 

XlL Station number 

SC Distance along the streamline 

sex Desired distance along the streamline. SCX Is proportional 

to the station number. 
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LC Station table Index 

lOOFC ' Indicator for primary stations (1) and all other 

stations (2) 

KCL K Index for the control streamline 

Values along orthogonal lines (K = l.NK) 

PHI2 Angle of the orthogonal with respect to horizontal; 

discrepancy In orthogonality. 

DSl Point movement along the streamline 

ZK Axial distance, abclssa 


RK 

WEZPT 

DSIC 


Radial distance, ordinate 

Double streamline Indicator (slip lines, wakes, stagnation 
streamlines). WEZPT(K) = 1 If K,K+1 is a double 
streamline. 

Correction to DSl to account for streamline curvature 
(negative [concave] curvature only ) 


Subroutine ; REFINE 

Arrays of data on streamlines (when new streamlines are being added) , or 
along orthogonal lines (when new stations are being added). 

Parameter Description - 


lA 

IB 

CR 

DELS 

DELVM 

LSTA 

MJ2 

SGX 

RAV 

ZAV 


Beginning of refinement 
End of refinement 

Criteria for refinement; Cj = Rg + Rg^^^R^ 
where Rg and Ry are defined in Section 12. 

Interval size of unrefined grid 

Velocity change for the unrefined Interval normalized by 
VMGl or VMG2 

Station Index 

Field point index (M) along streamline J2 
Grid size Interpolated from GZ&SGZ 
Average R of intervals to be refined 
Average Z of intervals to be refined 
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The error print Is best Illustrated by an example of the typical output 
after an error In the flow balance Coverlays (2,0), (2,1)). Referring to 
the example, the Initial output consists of the data In common /CFB/, Includ- 
ing a variable Identification key. This output Is followed by the Information 
In the streamline table commons /ALLCOM/, /CIDEX/, /IXORIG/, and, the STC 
tables BDYTAB, CONVTB, WAKETB. CADJ WF . and STATAB. (See Section 6 and 7 for 
definition of these commons and tables. T The field. tables are printed next 
and Include the following Items: 

Streamline number 
Field point Index 
Upstream field point 
Downstream field point 
ISTAG value 

Cumulative distance along streamline 
Ctmiulatlve distance along orthogonal line 
Axial coordinate, abclssa 
Radial or transverse coordinate, ordinate 
Flow angle (radians) 

Streamline curvature 
Velocity 

Coefficient B of matrix equation 
Right hand side of matrix equation 
Streamline adjustment 

The field table output Is followed by the /ERASE2/ data as described 
In the preceding section. 

9.2.2 Location and Explanation of Specific Error Comments 

As mentioned previously, Input or calculation errors will result In 
the termination of the problem. The locations (subroutine name) and descrip- 
tion of error conditions are given In the following table. Included are 
suggested corrective measures (where possible), description of comments 
preceding the table print, and a designation of whether the error Is FATAL 
or NON-FATAL. Several specific problems deserve further comment. When the 
Iterative solution of the matrix equations for streamline point movement 
Is In trouble (NSWP ->• MAXSWP) , a print of the solution matrix (DS2) and 
the maximum DS2 will start appearing in the printout. Subroutine lAD Is 
structured to detect this problem and provide the iteration history over a 
complete cycle as the solution proceeds from over-relaxation to under-relaxa- 
tion. In this case, the problem may usually be detected by inspection of 


J 

M 

MU 

MD 

I 

51 

52 
Z 
R 

PHIl 

CURV 

VM 

B 

RHS 

DS2 
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the solution matrix. Appropriate corrections would be to input a BHOBAS larger 
than the preset value or preferentially sweep the field in a single direction. 

The boundary layer data tables are built by subroutine BLTBBL at the 
end of the problem. If a boundary layer on a given surface is specified in 
the boundary layer input table and there is insufficient storage in the 
TABLES region to accommodate all the data, the following output will occur: 


TABLE SPACE EXHAUSTED-BOUNDARY LAYER DATA 
IlPPFTl 

FOR BOUNDARY "boundary name" NOT SAVED 

The above comment serves as a warning to the user to increase the 
available size of the TABLES area. Following the printing of this comment, 
the calculation switch in the input table is turned off. 

If subroutine BDYPTM detects a separated boundary layer in the course 
of interpolation for the displacement correction, the following comment 
with boundary name and separation location appears each time the bovmdary 
is accessed: , 

** WARNING** SEPARATED BOUNDARY LAYER, BOUNDARY = '*boundary name" 
SW = XXXXX 


Since displacement thicknesses downstream of the separation point are 
in error, the user is advised to terminate the problem at this stage. 

Several informative comments may appear in the printout during the course 
of the calculation. The specific diagnostics are as follows: 

• *** THE INPUT GRID REFINEMENT CRITERIA HAVE NOT BEEN SATISFIED 
appears after the Iteration history output if the specified grid 
refinement criteria have not been satisfied. 

• *** THE SOLUTION HAS NOT CONVERGED TO THE INPUT INTOLERANCE 
appears after the Iteration history output if the flow balance 
tolerance has not been satisfied upon completion of NINNER inner 
iterations. 

• ISTAG =3 POINT INSERTED L-E. OR CORNER AT (Z-loc. ) (R-loc. ) 

appears interspersed in the Iteration history output when the 
velocity at the stagnation point becomes less than half the velocity 
at the adjacent point above/below the discontinuity. At this stage, 
the stagnation point is converted to a "hard" discontinuity. 
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1 f 

o • u .1 t • •' U 

LA 1 







18 

8.096875 

EXT 

.969091 






19 

H. 062500 

EXT 

1.285959 






20 

8.125000 

EXT 

2.570908 






21 

8.250000 

EXT 

5.191817 






2? 

8.500000 

EXT 

10.283639 






23 

9.000000 

EXT 

20.567267 






29 

10.000000 

EXT 

91.139535 






25 

11.000000 

EXT 

61.701802 






26 

12.000000 

EXT 

82.269070 






27 

19.000000 

EXT 

123.903609 






28 

16.000000 

EXT 

169.538139 


- 



ALLCOM 






. 



1 

.800000 

19.695990 

518.686000 

1.000000 

585.080069 

-0 

1716.200 

1.900000 

9 

.800000 

19.695990 

518.688000 

22.901609 

585.080069 

0000 

1716.200 

oooo 

17 

0 

1.000000 

0 

..-0 

- • 

— 

-- 


CIDEX 









1 

300 

28 

283 

307 

1 

— 



IXORIG 







— 



1 

1 

20 







3 

21 

231 







5 

232 

269 




_ 


. . 



7 

270 

269 







9 

270 

269 

.. 






11 

270 

1173 







BOY T AB 









21 

CNTLN 

12 

0 

. . 

— 

* 

■ ■ 


2“ 

W2 

0 

0 

, 






■ . .. . 

27 

-30.000000 

0 

0 






30 

18.000000 

0 

0 


— ■ 


... . . 


33 

NACAl 

186 

6 






36 

W2 

EXT 

57 






39 

NACAl 

6 

57 

. . 





... . 


*2 

CLEX 

57 

177 






95 

18.000000 

8.900000 

3,207758 


_ ... 




. - , 

AB 

16.200000 

8.262330 

3,228091 






51 

19.900000 

8.091780 

3.291809 





._ 


59 

10.800000 

7.733290 

3.231861 






57 

8.100000 

7.529710 

3.199538 

_ 






60 

6.300000 

7.998560 

3.173916 






63 

9.500000 

7.913000 

3.159220 






66 

2.500000 

7.378000 

3.158959 






69 

.761000 

■ 7,398000 

3.127089 






72 

.539000 

7.361630 

3.025877 






75 

.383000 

7.387280 

2.913286 






78 

.198000 

7.951810 

2.679988 






81 

.115600 

7.501360 

2.512733 






89 

.083200 

7.527230 

2.919897 






87 

.059090 

7.556030 

2.298992 






90 

.037900 

7.575930 

2.203961 






93 

.017210 

7.610000 

2.023395 






96 

0 

7.682000 

1.570796 






99 

.001908 

7.696762 

1.316923 
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SECTION 10.0 


OPERATING PROCEDTJEES 


The STC program described herein may be nin on any Control Data 6A00/ 
6600 machine operating under SCOPE 3.0 or a higher level operating system. 

In general, operating procedures and control card set-ups will differ from 
site to site. The following comments on program modifications, deck set-ups, 
and operating Instructions are restricted to the program as Installed at the 
Langley Research Center. Minimal changes should be necessary for successful 
Installation at other CDC sites. 

10.1 GENERAL OPERATING PROCEDURES 

The version of STC available at the LRC computer site allows a maximum 
flow field data storage of 768 grid points, and requires a field length of 
115g K for execution from an absolute overlay file. 

A larger version of STC may be created by recompilation (FORTRAN IV RUN 
compiler) of the two block data subprograms which set the size of the STC 
table and field point storage (USECDG) and the size of the Infliience coeffi- 
cient arrays used in the solution for the streamline adjustment (USECDM) . 


1.) Block Data USECDG 

(overlay 0,0) 

Common 

Variable 

Tab le /Des cript ion 

/CDS2/ 

C12 

Field Table - Streamline adjustment 

/CRHS/ 

C16 

Field Table - Right hand side of matrix 
equation. 

/CHDATA/ 

C9 

STC Table Storage 

/CCURV/ 

Cll 

Field Table - Streamline Curvatures 

/CPHIl/ 

CIA 

Field Table - Flow Angle 

/CSl/ 

C17 

Field Table - Distance along streamlines 

/CS2/ 

C18 

Field Table - Distance along orthogonal 
lines 

/SLTAB/ 

C8 

Streamline Table - Dimension 3* maximum 
number of streamlines 

/CM/ 

C13 

Field Table - JMS (mesh point connection) 
array 

/CB/ 

CIO 

Field Table - Coefficient B in matrix 
equation 

/cz/ 

C20 

Field Table - Axial coordinate (X,Z) 

/CR/ 

CIS 

Field Table - Normal coordinate (Y,R) 

/CVM/ 

C19 

Field Table - Velocity 
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As Indicated in the preceding list, 11 variables comprise the flow field 
tables. These arrays should be dimensioned to the maximum manber of mesh 
points. The exception is /CDS2/, which must be set to a minimum of 900 
locations . 

Note also that /CZ/ and /CR/ must reside next to each other in memory, 
as these quantities are used to determine the maximum limit for the number 
of field points. For problems with a large quantity of boundary input, many 
boiaidary layer surfaces , and/or a large nuiober of anticipated orthogonal 
stations, it may be necessary to increase the size of the STC Table storage 
common /CHDATA/ beyond the 2200 locations currently allocated. The stream- 
line storage /SLTAB/ is set to accommodate 128 streamlines; this should be of 
adequate size for most problems. As in the case of /CZ/ and /CR/, the commons 
/CHDATA/ and /CEND/ must be located adjacent to one another in memory. These 
items are used to set the maximum limit for the STC table storage. 

2.) Block Data USECDM (Overlay 4,0) 


Common 

Variable 

Description 

/CA2/ 

C42 

Curvature influence coefficient 
A2 

/CA3/ 

C43 

Curvature influence coefficient 
A3 

/CA4/ 

C44 

Streamline flow difference influ- 
ence coefficient A4 

/CA5/ 

C45 

Curvature Influence coefficient A5 

/CA6/ 

C46 

Curvature influence coefficient A6 

/CA7/ 

C47 

Streamline flow difference Influence 
coefficient A7 

/CA8/ 

C48 

Streamline flow difference influence 
coefficient A8 


The dimension of the arrays in these commons should be the same as the 
dimension of the flow field arrays . 

10.2 RECOMMENDED PROCEDURES FOR PROGRAM MAINTENANCE AND MODIFICATION 

As mentioned previously, a source and absolute binary copy of the STC 
program are stored on the permanent disc file (data cell) at the LRC computer 
center. It is recommended that relocatable binary copies of each version be 
maintained on tape. As changes are made to the program, only those routines 
which differ need be recompiled. The relocatable binary tape may be updated 
using standard COPYL and COPYN functions. Subsequent execution from this 
file will result in the creation of a new absolute binary file which may 
replace the one currently on the data cell. The updated source decks may 
also be updated at this time. 
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The source copy of STC contains *DECK cards as the first card of each 
subroutine. In running at other .CDC sites, the source file may be used 
directly to Initialize an update tape. 

10.3 STANDARD EXECUTION DECK SET-UPS 

Normally, the absolute binary (overlay) program will be reloaded from 
the data cell for execution. The suggested deck set-up Is: 

V 

(Job Card) 

(User Card) 

FETCH (A3727, address on data cell, BINARY, STC) 

COPYBR (INPUT, TAPES) , 

STC. 

789 (EOR) 

[Input Data] 

6789 

As indicated previously, this deck set-up is applicable at the NASA Langley 
Research Center computer site. At other Installations, the absolute binary 
file may be reloaded from tape. For example, a typical deck set-up for use 
on the CDC Cybernet System Is : 

V 

(Job Card) 

LABEL (STCOVLY, R, VSN « UT6 11) 

COPYBF (STCOVLY, STC) 

REWIND (STC) 

COPYBF (INPUT, TAPE 5) 

STC. 

789 (EOR) 

[Input Data] 

67 gg (EOF) 

10.4 RESTART OF STC PROGRAM FROM DATA TAPE 

A partially completed STC problem may be restarted using the output data 
tape from the previous STC execution. The restart tape is obtained by specify- 
ing TAPOT as T on the general input header card and using a file card such as 
REQUEST or LABEL to assign TAPE2 to a given physical devlcie. Of course, the 
logical file TAPE2 must be assigned with "write permit" RING option. In the 
subsequent restart task, TAPIN must be specified as T on the general Input 
header card. Also, the Input file TAPEl must be defined using an appropriate 
file card (NORING) . The input for the restart case should consist of the 
Identification cards, the general input header card, and the first $A namelist. 
For example, to restart a job after 3 major grid refinements and run to 5 
grid refinements, the example input ml^t appear asi 
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NAME = J. Smith 

ADDRES = LRC 

IDENT = NASA 8 INLET 

2 4 14 24 

1 STC T F 

$A 

MAXIT =5, 

$ 

Intermediate output at a given MAXIT level may be obtained by nmning a restart 
case and not assigning TAPE2 and TAPEl. In this instance, these files are 
assigned to disc and switched at the end of each total set of input. This 
option is also useful for changing program tolerance at a given stage of the 
calculation. 

The preceding example could be nm as one case using the following input: 

V 

NAME = J. Smith 

ADDRES = LRC 

IDENT = NASA 8 1 INLET 

24 14 24 

1 STC F T 

$A 

MAXIT = 3 

• 

$ 

1 STC T F 

$A 

MAXIT = 5, 

$ 

Successive restart cases may be stacked in sequence. On all but the final one, 
use 1 STC. T T as the header card. 

As indicated in above, a partially completed STC problem may be restarted 
using the output data from a previous STC execution. When xised for a boundary 
layer run, the output tape may also contain the boundary layer data from the 
previous run. If TAPOT and TAPIN are specified as T and TAPE2 and/or TAPEl 
are not assigned to a tape file, the system will assign them to disc. Hence, 
consecutive botmdary layer restart cases may be run by simply setting TAPIN 
and TAPOT to T on all cases after the first. The general procedure for 
carrying out a boundary layer iteration woTild be to run a given STC problem 
to a specified refinement/or convergence level and then run successive restart 
cases at the same level to converge the combined inviscld - botmdary layer 
problem. 



SECTION 11.0 


HELPFUL USER HINTS 


11.1 GENERAL COMMENTS 


The following section Is a compilation of Information for the user. 

Many items have been covered elsewhere, but they will be repeated in this 
section for emphasis. The Intent is to, identify key input parameters so 
that computer solutions can be accomplished with minimum user problems or 
errors. 

The primary items of Importance are: 

1. Smooth input geometry 

2. Grid refinement definition 

3. Iteration tolerances 

4. How to get out of trouble when the solution doesn't converge. 

11.2 SMOOTH INPUT GEOMETRY 

The Streamtube Curvature method uses the boundary or body surface 
curvature as the geometry parameter which causes velocity and/or pressure 
gradients. Thus, the computer program is very sensitive to changes in cur- 
vature. Hence, an accurate definition of the boundary coordinates and local 
angles is imperative if surface curvatures are to be smooth and continuous. 

There are two methods to Insure smooth input geometry; Ci) use an analytic 
boundary or function with continuous first (angle) and second derivatives 
(curvature) or (2) smooth specified boundary coordinates By fitting a curve 
to the surface points. The definition of an analytic boundary is not always 
possible and, in some cases, the intersection of two analytic functions will 
have discontinuous second derivatives. (An ASME nozzle is an ellipse followed 
by a straight line . The discontinuous second derivative causes a pressure 
"Blip".! Hence, a curve fit technique may be necessary in most cases. 

The best curve fit to use is a piecewise smooth cubic since the curva- 
tures are evaluated in STC using the same type of curve fit. Hence, the 
definition of the coordinates and angles should come from some such geometry 
smoothing method. As part of the geometry definition, the curvatures should 
be checked for discontinuities or erratic behavior. 


11.3 GRID REFINEMENT 

The grid refinement criteria is explained in detail in Section 13. Some 
comments on use are appropriate with reference to Figure 31. 
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First of all, the purpose of the grid refinement criteria is to maximize 
the grid refinement in areas of interest and to minimize the flow field grid 
refinement in areas of smooth flow. By forcing flow field grid points to be 
used efficiently, computer table size and storage space is minimized and 
computer time is reduced. Also, the flow field grid can develop without 
over-constraining the streamlines. 

The best technique is to first Identify the region or regions in the 
flow field where a high degree of resolution is necessary. The geometric 
grid refinement criteria in this are^ is set at a value equal to the grid 
size desired. As shown in Figure 31, a rectangular flow field section on 
the inlet lip is defined with SGZ = SGR. If a subsidiary plot of SGR vs 
GR and SGZ versus GZ are set up, the input table is readily developed. In 
the radial direction, SGR is specified at the centerline, decreases linearly 
to the region of interest, increases linearly to the outer region value and 
is constant from there out. In the axial direction, a similar plot can be 
drawn. The coordinates of these two plots become the input tables for 
GZ/SGZ and GR/SGR. 

In some problems, the geometric grid refinement must be redefined if 
solution difficulties occur. A demand for excessive grid refinement in a 
local area can over-constrain the development of the flow field. Only experi- 
ence will allow the user to make the optimum choice of grid refinement 
criteria for different problems. 


11.4 ITERATION TOLERANCES 

The iteration tolerances, T0LINR and T0LES2, have been preset at 
recommended values for most problems. As explained in Section 9.1, these 
tolerances control the iteration logic for the grid refinement and the "flow 
balance" or inner loop solution. 

Usually, if the input geometry is properly defined and no errors are 
made in defining the flow properties (such as setting VARY = T in all the 
channels) , the solution proceeds through several grid refinements with no 
problems . 

When the grid refinement criteria and the iteration tolerances are 
incompatible, l.e., the tolerances are very small and the grid refinement 
has high resolution In an area of large streamline curvature, the solution 
will show instability in the definition of streamline adjustment (MAX-ES2 will 
not decrease continuously) . The remedy is to change the tolerances or the 
grid refinement or both. 

One good procedure is to keep the iteration tolerances relatively large 
during the flow field grid development and then restart the last grid refine- 
ment with smaller tolerances. Again, experience in use of the computer 
analysis will be necessary to provide guidance in how to do this. 
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11,5 TROUBLE SHOOTING 

The Streamtube Curvature Analysis has been set up as a user oriented 
computer program for the solution of relatively complex subsonic and transonic 
flow field problems. Many of the control variables have been preset to the 
most generally useful value. Hence, a majority of problems will be solved 
quickly and accurately if the Inputs (geometry and channel flow properties) 
have been defined, correctly. But, there will be instances when numerical 
difficulties will occur. 

Invariably, a large percentage of the numerical difficulties can be 
related to errors in defining the geometry or bad judgement in selecting the 
grid refinement or Iteration tolerances. However, when large supersonic 
bubbles are present in a transonic flow, convergence problems in the inner 
loop have been experienced. A dump, very similar to that shown in Section 
9.3, will result. The solution history gives the first diagnosis in that 
the level of grid refinement, the number of inner iterations during the 
last refinement, and the behavior of MAX-ES2 can be quickly reviewed. Often 
MAX-ES2 will start to converge and then show instability. By identifying 
the appropriate NINNER, the solution can be forced to the preceding level of 
grid refinement by rerunning the problem with the NINNER control inserted. 
Convergence difficulties for this type of problem may sometimes be overcome by 
establishing a refined grid at a lower Mach nxunber and slowly raising the 
Mach number to the desired level using restart cases. 



SECTION 12.0 


SAMPLE CASES - PROGRAM OUTPUT 


Three sample cases have been selected to illustrate the typical input/ 
output for the STC program. The first sample case is the analysis of an inlet 
lip at a free stream Mach number of 0.8 (NASA Inlet No. 8). This exanple 
case is’ identical to the one presented in the previous STC-SAB user manual, 
except the inlet surface has no boundary layer. The geometry is defined as a 
straight centerline for the lower boundary of the inlet flow (BUY CNTLN W2 
and UPPER = F) , the coordinates and angles for the inner surface of the inlet 
lip, the NACA Series 1 external cowl lip (CLEX) , and a far-field (FF) upper 
boimdary. Three (3) passes through the program are shown with MAXIT = 8. On 
the third pass, the tolerance (T0LES2) was reduced from 1. to .001 to obtain 
more accuracy in the invlscld solution. Also, on the final restart, the 
output indicator PRPRN was set to 1 to obtain a full printout of the flow 
field. 


The second case illustrates the bomdary layer on a two-dimensional 
circular arc airfoil in a wind tunnel at Mq = .663. As indicated in the out- 
put, boundary layer separation occurred at SW = 19.88. Upon restart, the 
previously described warning comment is printed each time the boundary is 
accessed to determine an orthogonal - boundary Intersection. 

The final test case Illustrates the revised flow adjustment procedure and 
alteration of the STC wake table to reflect boundary layer displacement effects 
at a trailing edge. The problem configuration consists of an outer nacelle 
portion (BODY) and an ideally expanded CD nozzle (NOZZLE, CNTLN) at a subsonic 
Mach number of 0.2. The nozzle lip trailing edge thickness (TTE) was 0.005 
and boundary layers were specified on both surfaces. The outer upper surface 
was taken as a far-field (FF) boundary. The solution was run to MAXIT = 8 
and a tolerance (TOLES 2) of 0.005 in two passes, the restart portion being 
used to provide displacement correction of the initial invlscld solution. 
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MATRIX SOLUTION PARAMETEPS- 

lADM = 0 1 » FOR STREAMLINE* ALTERNATING* AND ORTHOGONAL LINE RELAXATION) 

PH0RAS= .500 (ACCELERATION FACTOR* BASE LEVEL) 
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MATRIX SOLUTION PARAMETERS- 

lADM r 0 (=-I.0.I. rOR STREAMLINE. ALTERNATING. AND ORTHOGONAL LINE RELAXATION) 

PHOPAS= .500 (ACCELERATION EACTOR. BASE LEVEL) 

- RH0AMP= .500... (ACCELERATION EACTOR. AMPLITUDE OE VARIATION) 
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» * C A R 0 


INPUT** 


NAMF= DAVE FERGUSON 

ADDPES= . EVFNDALE. 

1DENT= LAHTI TEST CASE DEC. 1ST.1972 

. 1 STC F T... . - 

$A 

MACHO*. 6G3« 

TS0=537.736. 

_.PS0=1 0.9425. _ . 

AXI=F, 

-HAXIT = 5. 

RG=1716.2.VMG1=100..VMG2=100..SCON=198.6. 

SGR(1 )=1,0. 

PPPFN=-1, 

... S . . 

? ROY WALL FLOW 

...5A 

UPPFR=T.7R0NLY=F. 

. 8(1 )=-16. 0.9. 0.0.0. 

-9. 0.9. 0.0.0. 

— -8,0.9.0085. 1.022. .... 

-7,0.9.0373.2.327. 

_,-6,0.9.0905.3.703. 

“5 , 0 » . 1 f)4 0.4. 872. 

_-4, 0.9. 2535. 5. 502. 

-3,0.9.3503.5.273. 

0.9. 4300. 3. 744. __ _ _ 

-1 ,0.9.4812,2.139. ■■■ 

._0. 0.9, 5053, 0.887, 

1.0. 9.5072,-0.383. 

2.0. 9.4930,-1 .237. 

3.0. 9.4853,-1.857, 

.4.0.9,4303.-2,153. 

5.0. 9.3887,-2.880. 

_ 6. 084. 9,3344,-3.0591 

7.0. 9.2887,-2.725, 

_8. 0.9, 241 7.-2.474, 

9.0. 9.1991,-2.419, 

_.l0,0.9.1572,-2.372, 

11.0. 9.1171,-2.192, 

_ 12,0.9.0815,-1 .860, 

13.0. 9.0526,-1.447, 

14.0. 9.0309,-1.038, 

15.0. 9.0180,-0.685, 

_ 16,0 .9.0086.-0.405,, 

17.0. 9.0015,-0.182, 

_ 18,0.9,0,0.0,.... 

22.0. 9.0.0.0, 

.$ 

2 Rnv BUMP FLOW ' 

SA 

upper=f,zronly=f, ' 

„BL = T, . 

R(i>=-i6»o,o.6,o.b, 

. .■"6,064,0.0,0.0. 

”6»0.»0005,0. 844 , 

■5. 5, *0347, 6. 801, 

-5,0. ,1165, 11. 551. 



--..c-.37s5.ih.?:r/, 

- 3 . '=•. 5196 . 15 . 567 , 


-3 

-3 

-? 

-2 

-1 

-1 


. 364 ., . 556 «, 15 . 038 , 
. 0 «. 6492 , 13 . 4 . 05 , 
. 5 . . 7577 , 11 . 086 , 

. 0 • . 84.53 , 8.810, 

.5, .9129,6.594., 
.0..961 1 ,4..4.05^ — 


-0.5, .9900,2.210, 

- 0.0..9997.0.0-. 


0.5. .9900, -2. 211, 

1.0. , 96 1 1 , -4.. 406^ - 

1 .5. . 9129,-6.594, 

. -2.0.. 8453. -8.810, 

2.5. . 7577,-1 1 .086, 

.. 3.0, .6492,-13.405, 

3.364. .5568,-15.038, 
- 3.5, .5196,-15.567i 

4.0. .3755,-16.237, 

.4.5. .2350,-l4.R2L, 

5.0. .1165, -11. 550, 

_ 5.5, .0347,-6.797, 

6.0. .0005,-0.841, 


,-_6. 0 64 ,0. 0 , 0.8j 

22 . 0 , 0 . 0 , 0 . 0 , 

_s 

3 CHN TLOW 



VARY=r, 

„PG = 1716.32, 

TS0=537.726, 

._PSO= 1 0 .5A25j 

MACH0=.663, 

_1 

STC T 

_J>A 

MAXITsS, 

% 


T 



executing progm=stc 

TAP1N= r TftPOT= 


oLrof'-cor«‘ 0 »-'Ocr*-«i/iLnir»voruris£)«o^>j-r^r-“r*^Lr 4 'Ooo 
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lOENT* LAHTI TEST CASE DEC. 1ST, 1972 
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I0ENT= LAHTI TEST CASE DEC. 1ST. 1972 
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DSlnMP= .020 (STREAMWISE RT MOVEMENT DAMPING. =0 FOR NO DAMPING) 

DS10P1= .500 (ADDITIONAL STREAMWISE DAMPING ON FIRST PASS ONLY) 

HOOENS= 0 (REFINEMENT LEVEL TO WHICH CONSTANT DENSITY IS ASSUMED) 

RHOC = I. 000 RHOW = I. 000 RH0CSS= 1.000 RH0WSS» 1.000 (CORRECTION EO. DECEL. FACTORS) 



MATRIX SOLUTION PARAMETERS- 

laOH = 0 (=-1.0.1. FOR streamline. ALTERNATING. AND ORTHOGONAL LINE RELAXATION! 

RHOPaS= .SOO (acceleration factor, base LEVEL) 

PHOAMpr .SOO (ACCELERATION FACTOR. AMPLITUDE OF VARIATION) 
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LOWEP ROUMOiRY TO CHN=FLOW t STREAMLINE COORDINATE. X!2» .000 
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. A -R, N l..N_6 — SEPARATEO BL SOUNDARY= BOMP SW= EO. 118066 



243 


separated BL » BOUNOARY= BUMP SW= 20.118066 



separated 0L « BOUNDARYs bump SW= 20.118068 


rvj 

r> 
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separated BL ♦ ■ 8OUN0ARY= BUMP SW= 20.118068 



separated BL . ROUNOARV= RUMP SW= 20.118068 



245 


separated BL * eOUNOARY= BUMP SW= 20.118068 



separated BL ♦ BOUNOARY= BUMP SW= 20.118060 



246 


SEPARATED BL . ROUNOARY= BUMP SWa 20.11S06B 



SFPARATFO Bu * BOtJNOARYa BUMP SW« ?0. 1 18060 




247 


SEPARATED BL . BOUNOARY= BUMP SW= 20.118068 



separated BL t BOIINOARY= BUMP SW= 20.1 18068 
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IDEM= LAHTI TEST CASE DEC. 1ST. 197? 
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_0Sin.MP= .cao . .. (STREAMWISE RT movement damping. »0 for no DAMPING) 

nSI0Pl= .500 (ADDITIONAL STREAMWISE DAMPING ON FIRST PASS ONLY) 

...n00ENS= 0 .. (REFINEMENT LEVEL TO WHICH CONSTANT DENSITY IS ASSUMED) 

R^OC = 1.000 RHOW = 1.000 PH0CSS= 1.000 RHOWSS* 1.000 (CORRECTION EO. DECEL. FACTORS) 



-MATRIX SOLUTION PAPA“ETERS- 

= 0 (=-1.0.1. FOR STREAMLINE. ALTERNATING. AND ORTHOGONAL LINE RELAXATION) 

PH0HAS= .500 (ACCELERATION FACTOR. BASE LEVEL) 
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tOWER BOMNDA^'V PRESSURE -FORCE = - -.0153 ■ 

UPPER BOUNnARY PRESSURE FORCE = -.0000 

. SUB OF ABOVE A1.3518 

LEAVING momentum = A1.3671 

- EPPOP -- • -.0153 
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lOENT. REFLEXEO AFTER-BODY SUBSONIC NOZZLE 
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IOEnTs REFLExen after-body ^SUBSONIC NOZZLE 
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LOWER ROIINDARY pressure force =- -.0147 

UPPER ROIINOARY pressure FORCE = -.0000 

SUM OF ABOVE - 41.3525 

LEAVING MOMENTUM * 41.3h71 

error - » -.0147 
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SECTION 13.0 
PROGRAM INPUT SHEETS 
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STREAMTUBE CURVATURE PROGRAM 
WITH BOUNDARY LAYER 

Page of 

Overall Input Data STC/Sheet-1 

output tape? 

T or F 


Mach number, ambient pressure and temperature, fluid properties 
$A (1.) (1.) (1.) (1.4) 


MACH0= . TS0= 

, PS0» 

II 

. GAM= 

Hlghllgiht radius, maximum body radius 
RHL= , RM= , TTE= 

, body closure 

9 

tolerance 

» 

Axlsymmetrlc or planar? 
(T) or F 

AXI= 

spaclal grid refinement 
GR(1) = 

criteria, see 

> » 

notes 

9 

SGR(1)= 

» » 

9 

» 

NGR= 




GZ(1)« 

9 9 

9 

» 

SGZ(1) = 

9 9 

9 

9 


NGZ= , 

maximum Mach number Increment between grid points 


March 1975 


W 


Input tape? 
T or F 


1 STC 


streamwlse 

direction 

( 0 . 1 ) 


normal 

direction 


( 0 . 1 ) 


VMG1= _, VMG2= , 

maximum number of refinements, curvature damping. Incompressible start 

( 1 .) ( 0 ) 

MAXIT= , RH0C= , N0DENS= , 

fluid reference temperature, reference viscosity, Sutherland constant, 
vnlts conversion constant 

(518.7“ R) (lO'^lbm/ln.sec) (198.6“ R) 

TREF= , MUREF= SC0N=_ , 



STREAMTUBE CURVATURE PROGRAM 


March 1975 Page of 

STC/sheet-2 

Boundary Coordinates 


boimdary 

name 

channel 

name 



w w 

V 



2 BDY 




upper botmdary? 
T or F 

angle Input? 
T-no, F-yes 

boundary layer? 
T-yes , (F-no) 

equiv. flat plate distance 
to boundary layer origin 




(0.) 

$A ^ UPPER= 

, ZRONLY= 

, BL= 

. CAPX1= 

Z 

R 

ANGD 




$ 
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STBEAMTUBE CURVATTOE PROGRAM 

March 1975 Page of 

Channel Flow Properties 

channel 

name 

3 CHN 

$A . 

ratio of flow rate 

specific gas may be 

heats constant adjusted? 

(1.4) (1.0) (T) or F 

GAM= , RG= , VARY= , 

stagnation properties, see notes 3 and 4 

total temp • total pressure 

TT0=_ , PT0= , 

Mach no. static temp static pressure 

MACHO , ISO PSO , 

flow area normalized by Ag^ 

AO= , 

$ 
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General Instruction and Notes for Sheet-1 of the STC Input Forms 

1. The STC Program computes the subsonic and transonic field of inviscid 
flow past (and within) arbitrarily shaped planar and axisymmetric bodies 
Inlet and exhaust nozzle flows wherein there may exist jet streams with 
differing energies are typical applications. An optional boundary layer 
analysis (SAB) is Included to evaluate friction losses and displacement 
of the inviscid flow. Note - on Sheet-1 and the following input sheets, 
the values in parentheses are used if other values are not input . 

2. The total flow is composed of one or more streams, the properties of 
which are to be listed on Sheet-3 (except as noted below) . Each stream 
occupies a "channel" which la identified by a one to six character 
alphanumeric word. Each channel must be boimded, at least in part, by 
an "upper boimdary" and a "lower boundary." Each boundary is also given 
an Identifying one to six character name and the coordinates are listed 
on Sheet-2. The following sketch Illustrates the naming of channels 
(CHN) and boundaries (BDY) . 


.z: — ^ y y X y ^ ^ y y y 

I L bdy=ff 

CHN=EXT 



An external flow channel must be named EXT, and the recommended name for 
the inlet flow capture flow channel is INT. The far-field Interface 
boundary must be named FF. Up to two arbitrary pressure or free 
boundaries are allowed and must be named PRESl, PRES2, FREEl, FREE2 (see 
Note 12) . With these exceptions, the selection of channel and boundary 
names is arbitrary. The special channel names EXT and INT cause extra 
streamlines tp be placed in the first refined grid. The boundary name 
FF indicates that the boimdary condition on FF is to be obtained from 
an analytic far-field solution. 

Program limitations preclude configurations with non-staggered leading/ 
trailing edges (see Note ^) . If a multistream case is to be analyzed, it 
is preferred that the low total pressure stream be the outer (or upper) 
channel. There are no si>ecific limits on the number of channels or the 
number of botmdarles. However, the total amoimt of data which may be 
input is limited by the available memory. 
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3. The solution method consists of constructing a grid of streamlines and 
orthogonal lines. Starting with two streamlines per channel (one for 
each boundary) and an orthogonal passing through the first and leist 
point of each boundary, the grid Is automatically refined by dividing 
the grid Intervals In half and In half again as required. The ninnerlcal 
resolution, the solution accuracy and the computer execution cost are 
all directly related to the extent of grid refinement. The Input vari- 
able MAXIT determines the maximum number of refinements . Providing 
this limit Is not exceeded, the grid will be refined locally, as re- 
quired, until the spacing of orthogonals and streamlines is less than 
the value determined from the SGR and SGZ tables and the Mach number 
difference between any two points on a streamline or an orthogonal line 
Is less than VMGl and YMG2, respectively. Grid size values versus raditis 
(or y-ordinate) are to be tabxilated after SGR and GR, respectively. NGR 
is the number of entries in each list. Grid size versus the axial 
coordinate Is to be tabulated after SGZ and GZ, respectively. NGZ is 
the number of GZ values. If dimensional values of RG, TSO and PSO are 
input (See Note 6) , then VMGl and VMG2 must have units of velocity rather 
than Mach number. See supplemental notes for additional details. 

A partially refined grid may be saved on tape by specifying a T in 
column 24 of the first card, or read from a previously created tape by 
specifying a T In column 14. 

If TAPEl and/or TAPE2 are not assigned via a REQUEST card, they are 
assigned to disc. This allows the user to obtain output for a given 
refinement level and provides the option of changing input parameters 
on the restart. For the restart case, specify a T in colinnn 14 of the 
first data card and include In the $A list only those Input quantities 
(viz; MAXIT) which differ from those originally input. 

4. In the initial calculation grid, an orthogonal line will pass through 
each leading and trailing edge point and through each sharp comer point 
(l.e., a point on the boundary with an angle discontinuity). It is not 
possible to analyze a configuration in which two or more of these points 
are approximately opposite to each other. For example, if a configuration 
contains more than one leading edge, the edges must be staggered relative 
to the streamwise direction. 

5. A free stream Mach number is specified by supplying a value of MACH0. 

6. Perfect gas assumptions are employed and the levels of ambient pressure 
and temperature may be dimensionless (TS0=PS0=RG=1) or dimensional. 
Dimensional values in a consistent set of units as described in the STC 
User Manual must be supplied if boundary layer calculations are requested. 
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7. A reference (or highlight) area is calculated from the input value of 
R ht. as follows : 


axlsymmetric : 
planar: 


^iL “ ®HL 

“ ®HL 


This reference area (or Ay in the planar case) is used for defining the 
mass flow for each channel. See STC/Sheet-3 note 5. 

8. Computed pressure drag forces are normalized by the (maximum) body area 
where 

2 

axlsymmetric: A = ir R 

mm 

planar: A = Ay = R 

/ m in in 

9. Finite trailing edge thickness is permitted; the maximum thickness, or 
body closure tolerance, is to be siq)plled after TTE. 

10. In some cases, it may be necessary to Increase the curvature damping RHOC. 

RHOC is preset to unity; it should be input as 2 or more if an oscillation 

of MAX-ES2 is observed in the Iteration History Printout. Studies are 
needed to provide better guidelines for the selection of , the RHOC 
parameter. 

11. NODENS is the number of grid refinements for which the streamline posi- 
tions are foimd by using a constant density (based on the total tempera- 
ture and total pressure) . This option eliminates the possibility of 
velocities greater than sonic and is recommended to ensure reliable 
starting. Example: 

N0DENS=2 , 

The default value, N0DEN=0, causes the program to use the Incompressible 
equation for the zeroth refinement only. For jet plume analysis, where 
adjacent streams have large differences in total pressure, NODENS =-l is 
recommended. 


12. Three special boundary options are provided: far-fleld, arbitrary 

pressure, and free. The far-fleld option may only be used for an 
external flow upper boundary. This will cause the numerical solution 
to be notched to a small perturbation analytical solution in the region 
from the "far-fleld" boundary to infinity. This boundary option is 
invoked by specifying the boundary nan® FF (See STC/Input Sheet 2. - 
Note 1) . 
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An arbitrary static pressure may be specified along one or two bomdariesi 
The input is : 


NZP = 
ZP(1)= . 
PPS(1)= 
NZPl = 


> 

,> ^ > > » » > 9 9 9 

,9 9 9 9 9 9 9 9 9 

, PSPISV = , 


The allowable pressure boimdary names are PRESl and PRES2 (See STC/Input sheet 
2. - Note 1) . 

PSP and ZP are tables of pressure versus aiHal distance. If there are 
two pressure boundaries, the values of Z and PS for the second bomdairy are 
to be included in the same ZP and PSP arrays immediately after the values for 
the first boundary. Also NZPl must be defined as the number of values which 
apply to the first boimdary. For each pressure boundary there must be at 
least two ZP, PSP data values and ZP must span the region in which the boundary 
pressure is to be controlled. The total number, NZP, cannot exceed 10. 


If PPS input values are velocity rather than static pressure, then set 
PSPISV=1. For a jet bounded by two PRES or FREE boundaries, also specify the 
downstream angle: ACF(2)= , NBCIN(2)=1, (Refer to Note 13). 

A free boundary is a constant pressure boimdary which is downstream of 
a fixed boundary and the pressure level is the same as the pressure at the 
last point on the fixed boundary. Thus to use the free boundary option the 
contour must be input as two (or more) boundaries (see Note 2, Sheet-2). 

The allowable free boimdary names are FREEl and/or FREE2. 

The free boundary option does not converge easily since the change in 
boundary pressure from one Iterate to the next is not included in the stream- 
line position correction equation. 


For all of these options the approximate boundary shape must be input in 
the usual manner using Sheet-2. This is required to set up the initial solu- 
tion grid: the boundary will then "float" in order to satisfy the specified 

boundary conditions. For the PRES or FREE options, it is recommended that 
the grid be developed by using fixed walls for the first few refinements. 

This may be enforced by setting N0DENS=3, say, producing an incompressible 
grid and delaying the pressure boundary option to the fourth refinement. 


13. Boundary conditions on the streamlines at the upstream and downstream 
boundaries may be input as follows: 


Upstream: NBCIN(1)= , ACF(1) = 

Downstream: NBCIN(2)= , ACF(2) = 


NBCIN=1 if ACF is angle in degrees 
=2 if ACF is curvature 
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NBCIN and ACF are preset to 2 and 0, respectively. This yields uniform 
flow with no cross stream pressure gradient. The location of the 
upstream and downstream boundaries are determined by where the body 
surface, centerline and/or fair-field contours end. Since b.c. data are 
usually not known precisely, these contours should be extended (about 
one-half channel width) beyond the region of interest. 

14. The fraction of unblocked circumference or lamina thickness (in the 
planar case) may be input in a rectangular array using the following 
format : 

NTHKX = , NTHKY = , 

THKX(I) = , , _, . , , , , 

THKY(1) = , , , - , ' . , , 

THIK2D(1)= , , , , , , , 


> ♦ j ♦ » » » » » 

etc. 

THKX are the axial positions of the tabulated lamina thickness in the 
THIK2D table; list the y-variations first, i.e. , 

([THIK2D(J,D ,J=1,NTHKY], 1=1, NTHKX) 

Outside the range of THKX or THKY, the end values in the THIK2D table 
will be used. The THIK2D table is not extrapolated. Table size limits 
are as follows : 

2j<NTHKXj<25 1j<NTHKYj<25 NTHKX*NTHKY_<250 

15. At each level of grid refinement, a sequence of inner iterations is per- 
formed to converge the non-linear streamline position equations. Each of 
these inner iterations includes the evaluation of the streamline position 
error, ES2, and then a matrix solution of the linear correction equation 
for the streamline adjustment, DS2. The maximum nxnnber of inner iterations 
may be input as follows ; 

NINNER(1)= , , , , , , , 


Likewise the fraction of computed DS2 actually used when moving the stream- 
lines is : 

CNVF(1)= , , , , , , , , , 


In these tables, the index (which may go as high as MAXIT 16) is the 
refinement level. All values of NINNER are preset to 10 and all values 
of CNVP are preset to unity. NINNER and CNVF input values are seldom 
needed. 
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16. During the inner Iterations, the grid points are moved in the streamh^lse 

direction to obtain an orthogonal pattern. The following streanwise 

danq>lng factors are available but rarely needed: 

(0) (.5) 

DS1DMP= , DS1DP1= , 

DSIDPI is an additional factor on DSIDMP for the first inner iteration. 

Use values between zero (for no damping) and unity (for no movement) . 

17. Built in tolerances which may be changed by user input are as follows: 

T0LES2 .001 final solution tolerance on the streamline position error, 
MAXES2, as determined by the flow balance routine. The 
input value of T0LES2 is first multiplied by the character- 
istic length, defined as the average of the SGR and SGZ 
values . 

TOLWF .001 toleance on the fractional flow adjustment needed to meet 
the trailing edge pressure closure condition. 

TOLINR .05 intermediate solution tolerance on streamline position 

error, MAXES2. This tolerance is first multiplied by the 
current grid size, is printed as LIM-ES2 and is satisfied 
(or NINNER is exceeded) before proceeding to the trailing 
edge pressure closure iteration or to additional grid 
refinements. 

TOLRL .001 matrix solution tolerance on the sweep- to-sweep change in 
DS2 divided by maximum DS2. 

18. Miscellaneous input items: 

PRPRN = c 

= (0) for a, b & c 
= -1 for b & c 

code: a - flow properties at all grid points printed 

b - flow properties along boundaries printed. 

c- flow properties along approach and trailing stream- 
lines printed. 

CHOTST = T (preset value) 

= F to delete the check for choking when iterating for trail- 
ing edge pressure closure, set to False only for low 
speed flows to reduce execution time 
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19. Boundary layer calculations are normally specified using the input given 
on STC/Input Sheet-2. In some cases it is desirable to introduce bound- 
ary layers into a fully converged STC invlscld solution. This may be 
conveniently done on restart by specifying INPBL as the total number of 
surfaces for which a boundary layer calculation is desired. The boundary 
names and initial X's (CAPXl) are given on fixed format cards immediately 
following the $A NAME LIST. The format of these cards is as follows: 

Column 2-7 Boundary name (1-6 alphanumeric characters) 

Column 12-21 Equivalent flat plate distance from boundary layer 
origin to the first calculated boundary layer point 
(F10.6 format). 

Example: 

V 

$A 

5 Bestart Input 

INPBL=2, 

$ 

2 12 

BDYl Oq 

BDY2 1.632 

20. For boiBidary layer cases, supply a reference temperature, reference 
viscosity and a Sutherland constant if different from air values. The 
following sets of units may be used: 
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Units 


Parameter 

Dimensionless (STC) 

English (in.) 

English (ft .) 

MRS 

L 

any 

in. 

ft. 

m 

PSO, PTO 

*atm 

psla 

psfa 

~ N /m^ 

TSO, TTO 

*atm 

“R 

“R 

Or 

TREF 

- 

*R 

“R 

“R 

MUREF 

- 

Ibm/ln.sec 

Ibm/ft.sec 

kg/m. sec 

SCON 

- 

®R 

“R 

“R 

RG 

1 

ft2/sec^“ R 

ft^/sec2® R 

J/kg °R 

CG 

- 

ft-lbm/lbfsec^ 

ft-lbm/lbfsec2 

(unity) 

VMG1,VMG2 

** 

ft /sec 

ft/sec 

m/sec 


* 

atm - Normalized by ambient conditions 
** 

- Dimensionless (values are approximately equal to a Mach number 
difference) 


Notes for Sheet-2 of the STC Input Forms 

1. Use one of these sheets for each boundary. Supply a one to six character 

name to identify the boundary in column 14 of the first card. Special 
boundary names are: /FF - Far-field/, /PRESl, PFES2 - Arbitrary pressure/, 

/FREEl, FREE2 - Free/. (See Note 12 - STC/Input Sheet 1) . Also Indicate 
the name of the channel to which the boxmdary is adjacent in column 24. 

On the second card indicate whether the boxmdary is above (UPPER = T) or 
below (UPPER = F) the channel. 

2. The upper or lower "contour" which bounds a given stream may be composed 
of several "boundaries." In this case, an input sheet must be completed 
for each boundary; the last point of the first boundary must have the 
same coordinates as the first point of the second boundary, and so forth. 
This option is useful when considering variable geometry configurations 
such as flaps or movable nozzle parts. The movable part may be trans- 
lated and rotated, as indicated by Note 8, while the fixed part is held 
stationary. 

3. List values of Z (or X) , R (or Y) and the surface angle in degrees at 
discrete points along this bomdary contour after the symbol "B(l)=". 

Points at sharp comers must be listed twice, one time for each angle 
which exists at that point. In each interval, the STC Program fits 

a locally rotated cubic polynomial. The input points must be smooth 
and consistent with the specified angles. 
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All points are to be listed In the streanwlse direction. For an Inlet 
lip, the points are listed by starting at the hl^ll^t point and then 
proceeding around the nose to the trailing edge or downstream boundary. 

The Internal and external surfaces are listed separately under different 
boundary names. However, the coordinates of the first point must be the, 
same with ANGD equal to +90® for the external surface and -90® for the 
Internal surface. 

It Is recommended that the boundary coordinates and angles be obtained 
from an analytic definition of the contour, and that around the nose, 
angle variations between points be 20® or less. 

4. Pressure and Mach number distribution data will be printed at each 
orthogonal Intersection with the boundary, and not at each Input boundary 
point. Orthogonal stations, however, will be placed at any repeated 
point In the boundary table. List the san£ points twice If It Is desired 
to have an orthogonal placed In that position. (This option is modified 
slightly when ZRONLY = T.) Orthogonal stations are always placed at the 
beginning and end of a boundary and at a jimcture point between boimdaries 
along the same contour. 

5. If the coordinates but not the angles are known, the third column in the 
B-input array may be omitted. In this case specify ZRONLY = T and list 
the coordinates twice at any point where a curvature jump or an angle 
jump exists. The double points will later be deleted if the angle 
discontinuity Is less than 0.01 degrees. These double points are 
removed because extra calculation stations are usually not desired at 
such points. However, the double point angle tolerance, DBLPTS, preset 
as 0.01, may be input as zero if such double points are to be retained. 

6. With either input option, care must be taken to specify the coordinates 
with precision. The round off or reading error of the coordinate data 
should be less than As2/(10*L), where AS Is the local distance between 
points and L Is some characteristic length, say the length of the cowl. 
Conversely, the spacing between points should not be less than (10 6L)l/2 
where 6 is the relative accuracy of the coordinate data. The tabulated 
output curvatures may be consulted to verify the smoothness of the input 
data. 

7. NACA Series 1 Cowl coordinates are stored internally. With the ZRONLY = T 
option they may be selected by listing: 

B(l) = 991, 1, 

X]^, Y^, Series 1 Segment 

^ 2 . ^ 2 , 

X2. Y2 

_, , Cowl Aft Segment 
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where Xi, Yi are the hl^llght coordinates and X 2 , Y 2 is the position 
of the maxiimim diameter at the end of the Series 1 contour segment. 

8. The input coordinates of a botmdary may be adjusted by siipplying the 
following input quantities not shown on the front of this sheet : 


ROTATE 

ZPIVOT 

RPIVOT 

SCALE 


ZTRANS 

RTEANS 


angular rotation in degrees 
pivot coordinates 

multiplicative constant to be applied to the coordinate 
data 

translation increment in the axial direction 
translation on increment in the radial /vertical direction 


The order of transformation is rotation, scaling and translation. Hence, 
the pivot coordinates are in the same coordinate frame as the input data 
and the translation Increments are in the rotated coordinate frame after 
scaling. It is only necessary to input data for the transformation 
operations to be executed. 

9. The normal option assumes no boundary layer (BL ® F) . If a boundary 
layer calculation is desired, input BL = T. Also, supply an initial 
"equivalent" fat plate distance for the boundary layer origin to the 
first calculation station if different from 0. (stagnation point) 


r M , 

^ ^ 1+.2M^^ 

a = 1.25 Rex^ == 10^ 

1.2 Rex^ = 10^ 


CAPXl = 


±1 

Vl' 


/L^as 


orig 


Notes for Sheet-3 of the STC Input Forms 

1. Use one of these sheets for each channel to supply entrance flow proper- 
ties. (See exception under Note 5.) 

2. Of the input items shown on the face of the input sheet, use only those 
which are required for the selected options. 
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3. The total pressure and total temperature may be Input by either of the 
following two procedures: 

a. Specify TTO and PTO If the stagnation properties are known. These 
values may be normalized by the free stream ambient temperature and 
pressure. 

b. Specify MACHO, TSO and PSO If the static properties and Mach number 
are known. Again TSO and PSO may be normalized by the free stream 
ambient values. If only MACHO Is supplied (TSO and PSO are omitted), 
the TSO and PSO values from STC/Sheet-1 will be used. 

If neither of the above Is Input, free stream values as supplied on 
Sheet-1 are used for MACHO, PTO and TTO. 

4. If the gas constant, RG, Is different from the value supplied on STC/ 
Sheet-1, supply the value which applies to this channel. RG, TSO, TTO, 
PSO and PTO may all be given as dimensionless (normalized by free 
stream ambient) or dimensional using a consistent set of units. 

5. Input a value AO for the determination of the channel flow rate. AO is 
an area fraction normalized by A m. as defined under Note 7 of Sheet-1; 
the (dimensional) channel flow area Is then the product AO*AhL' Th® 
flow rate for the channel Is conq>uted by using one-dlmenslonal relations 
from the total properties (as determined under Notes 3 and 4), the 
supplied toch nu^er, MACHO, and the flow area. For Internal Inlet 
channels, specify R^l ^ the hl^llgiht radius and AO as the mass flow 
ratio. 

6. If for any channel the Input data on this sheet Is not supplied, the 
reference properties on STC/Sheet-1 will be employed and the frontal 
area calculated at the entrance station will be used. This option Is 
suggested for an external stream. 

7. Although approximate flow rates must always be supplied according to 
Note 5, the program will adjxist channel flow rates as required to meet 
the zero pressure loading conditions at a trailing edge or to meet a 
maxlmimi (choked) flow rate. The number of channels which require flow 
rate adjustment Is equal to the niimber of trailing edges. If the flow 
rate Is not to be varied for this channel, specify VARY = F. 
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Grid Refinement Criteria - Supplemental Notes 


One of the Important controls that the user should exercise In operating 
the STC program is the Grid Refinement Specification. In simpliest terms 
the user is required to input the following three quantities: 

SGR = , Nominal distance between streamlines and between 

orthogonal lines 

VMGl = , Nominal Mach number (velocity) difference between 

points along a streamline 

VMG2 = , Nominal Mach niunber (velocity) difference between 

points along an orthogonal line 

Note, if dimensional values of the perfect gas constant and tenq)erature 
are input then VMGl and VMG2 are velocity differences; otherwise VMGl and 
VMG2 are to be input as a Mach number difference. 


To Illustrate the function of these input items , consider the question 
of whether or not a new orthogonal line should be inserted between existing 
orthogonals as shown in Figure 30. To answer this question, the program 
computes, say, between the (i) «md (1+1) orthogonals, the following quantities 


(SG = SGR) (la) 


at each streamline: 
ASl 


R 




1 


SG 


V V 

\ l - i»jl 

j VMGl 


j = 


1,2,3, 


(lb) 


Figuratively speaking, these values are then plotted on Figure 31. 

If all the Rsj ,Rvj points fall outside of Region (1) a new orthogonal 
will not be inserted into the 1,1+1 interval. But if one or more points do 
fall in Region (1) , a new orthogonal is required. The new line will be 
extended to all streamlines • for which the Rsj ,Rvj values fall in Region (2), 
or to a minimum of five streamlines. 


Thus, according to the criteria Illustrated in Figure 31, the maximum 
distance between any two points along a streamline is SG (=SGR) . If a 
velocity gradient exists along the streamline, the maximum distance between 
points is reduced. When the spacing ratio Rs becomes very small (Rs < 0.24) 
relatively large velocity differences between points is allowed, a feature 
which is intended to eliminate excessive refinement near stagnation points. 
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Figure 31. Grid Refinement Criteria. 



A slnillar criteria is used for determining when additional streamlines 
are required. The only variation is that Equations (la) and (lb) are replaced 
by Equations (2a) and (2b); 



AS2i 

SG 


(SG = SGZ) 


(2 a) 


V V i 

\ - I - 1.J| 

1 VMG2 


1 = 1.2,3, 


(2b) 


Streamlines and orthogonal lines are inserted until the above criteria 
are satisfied (or until the number of refinements equals MAXIT) . The internal 
values of VMGl and VMG2 (0.1) are appropriate for most configurations, and 
unless these values are to be altered the VMGl and VMG2 input may be omitted. 
However, the user must always supply a valxie for SGR. In selecting this 
value care must be taken without sacrificing computing efficiency with an 
excessive number of grid points. 


As a guide to the selection of SGR it is suggested that SGR be set to 
about twice the expected grid spacing. Because new grid intervals are always 
obtained by halving existing grid intervals, the. average grid size will be 
less than that required to satisfy the criteria given in Figure 31. 


Spatial Variation of Grid Size 

For external flow problems, and some intemal problems as well, it is 
necessary to input a spatial grid size variation. Hie additional input 
required is ; 

GR(1) = , _, , , > 

SGR(1)= , , , , » 

GZ(1) = , , , , , , 

SGZ(1)= , , , , , , 

GR and SGR is a table of grid size (SGR) with radius (GR); likewise GZ 
and SGZ is a table of grid size (GZ) with axial position (SGZ) . The value 
of SG in Equations (la) and (2a) at any point, z,r, is then taken as the 
maximum of the SG's found by linear interpolation in the two tables of radial 
and axial variations. 

Special rules which apply to these input items are as follows ; 

1) If a radial variation (but not axial variation) of grid size is 
desired, the user must input pairs of values in the SGR and GR 
lists. GR must be in ascending order. If grid size information 
is required beyond the limits of the input table, the closest 
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values will be used. Also if the grid size Is to be constant 
with radius no GR Input Is required and only one value of SGR 
need be sixpplled. 

2) If an axial variation (but no radial variation) of grid size is 
desired, the user must supply pairs of valties in the SGZ and GZ 
lists . Also he must Input one val\ie of SGR which Is the minimum 
of the SGZ values. 

3) If both an axial and radial variation of grid size is desired, he 
should supply appropriate values in both tables. Remenber that 
for any point, r,z, the r-coordlnate will be used to obtain a grid 
size from the GR.SGR table and the z-coordinate will be used to 
interpolate in the GZ,SGZ table. The maximum of the two Interpo- 
lated values will then be used in Equations (la) and (2a) Generally, 
the SGR and SGZ values are set to small values close to a body, or 
region of Interest , and Increased to larger values as one moves 

away from the body. 


Input Control of "Partial” Grid Lines 

As noted previously, new orthogonal lines and streamlines are not 
necessarily extended to the field boundary. Instead, the length of these 
lines is determined by the size of Region (2) in Figure 31. Notice the 
Region (2) is Increased, and the length of the new grid lines is likewise 
Increased, if CRX, the constant in the equation defining the lower boundary 
of the region, is decreased. If CRX = 0, then all new grid lines will span 
the field; that is, full grid lines rather than "partial" grid lines will be 
Inserted. 

Values of CRX are preset, but they may be altered by input cards as 
follows : 


Input Name 

Preset 

Used for Extending 

(a) 

(b) 

Value 


CRX(l) 

CRXSL 

.375 

New Streamlines 

CRX(2) 

CRXOL 

.375 

New orthogonal lines across 
a subsonic region 

CRX(3) 

CRXSS 

.125 

New orthogonal lines across 
a supersonic or mixed flow 
region 

CRX(4) 

CRXE 

0 

New orthogonal lines which 
cross a sonic line 

CRX(5) 

CRXC 

0 

New orthogonal lines which 
cross a supersonic to sub- 
sonic compression line 

If two or more 

of the 

above values apply. 

the smaller value will be 


used. 
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Recommended Spatial Grid Variation for an Inlet Configuration 

Typical input for the spatial variation of grid size for an inlet 
configuration is given in Figure 32. In this case all dimensions are normalized 
by the highlight radltis. The quoted values of SGR, GR, SGZ, GZ should be 
scaled by the hi^ligjit radius for a given problem. Notice that the conbina- 
tlon of the axial and radial variations yields contours of grid size which 
are rougjily centered about the leading edge. 
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APPENDIX A 


LOWER BOUNDARY CORRECTION EQUATION FOR VELOCITY B.C. 


Continuity form of correction equation for first streamtube: 
(See Equation 94 of Reference 1) 


(pV)i^3 [6A2 - 6A^] + [A^2 - Aoi3 fi(PV)i .5 




( 1 ) 


We need now an expression for 6(pV)jL 5‘ From Equation 105 of Reference 
1 we have 


fi(pV)i^5 = ^1.5 Pi.5<5Vj_^5 


( 2 ) 


From Equation 95 of Reference 1 we have 


«1.5 “ «''l « 'I'l.S - V 


(3) 


SO 




(4) 


X 3V 

To evalxiate (Vi ^5 - use the cross stream momentum equation, ^ ~ 

Again following Equation 99 we can write 


^1.5 " ^1 “ “ 2 ^1.5 ^“2 " “1^ ^1.5 


(5) 


An approximate variational form of the above Is: 


'<^1.5 - V = 2 \.5 <“2 - “l> 


C6) 
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( 7 ) 


Here, we assume and that 


Vj^jCnz - ni) C^_5 


Equation [6] substituted into [4] gives 


- i ''i. 5‘”2 - ”l) 


( 8 ) 


And substitution of [8] into [1] gives 


3 

6A2 - 6A^ + [Aq2 ~ Aqj^] [ 6 V^ - I 3<n2 - n^ dC^J ( 9 ) 

X* 3 




Rearranging 


[6A2 - 6 A ^3 + 2 (^02 “ ®1.5 ^”2 " (10) 

6V^ 

~ ^o\ ~ ^o\ ~ ^^02 " '^01^ ^1.5 V . 
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